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ABSTRACT 


The directional properties of cosmic ray telescopes are analyzed with regard to both geo- 
magnetic deflection of the primary particles and to the spread in the atmosphere of the second- 
aries. 

The concepts “image of a cosmic ray source” and “‘polar diagram of a telescope outside the 
earth’s magnetic field’? are introduced. Image of a cosmic ray source is defined as the intensity 
variation recorded by a telescope during one day. The polar diagram represents the sensitivity 
of a telescope for different angular positions of a cosmic ray source element on the celestial sphere. 

A few examples of calculating polar diagrams are given. The resolving power of cosmic ray 
telescopes is discussed. 


INTRODUCTION 


The purpose of measuring the extra-terrestrial intensity variation of the cosmic 
radiation is to find the positions of cosmic ray sources in space as well as to study the 
influence of cosmical electric and magnetic fields on the particles and thus explore 
the electromagnetic state within the planetary system. All this can be looked upon 
as the problem of finding the shape and strength of a cosmic ray source located far 
away from the earth. 

The measurements are usually carried out with Geiger-Miiller tube telescopes, 
ionization chambers and neutron recording instruments. Telescopes measure radia- 
tion coming from directions within a rather narrow cone and have a higher resolving 
power than ionisation chambers and neutron recording instruments, which in prin- 
ciple are omnidirectional. Owing to atmospheric absorption, however, particles 
arriving in the zenith direction contribute more to the counting rate than particles 
from other directions, a fact which gives the omnidirectional instruments properties 
similar to those of broad angle telescopes pointing in the zenith direction. 

If observations are carried out using only one telescope with a well defined measur- 
ing direction the information derived when the earth rotates is limited to a narrow 
annular region centred around the earth’s axis. In order to get more information 
about the cosmic ray intensity distribution in space we must use several telescopes 
at different latitudes with each telescope scanning a certain region around the earth’s 
axis. This can also be achieved with telescopes pointing in different zenith- and azi- 
muth directions at the same place on the earth. The most suitable locations and 
directions of measurement depend on the deflection of the primary cosmic ray particles 
in the geomagnetic field. 
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Presentation of the problem and definition of cosmic ray image 


In general the cosmic ray source viewed from the earth extends over a certain region 
of the celestial sphere and the particle intensity varies over this region. Suppose that 
the angular extension of the source and the distribution of the particle intensity on the 
celestial sphere as well as the differential momentum spectrum of the particles are 
given. 

We want to know the intensity variation recorded by a cosmic ray telescope at 
a given place on the earth pointing in a given direction when the earth rotates and 
the telescope thus scans a certain region of the source as in fig. 1. 

The intensity variation recorded by the telescope during one day is called the 
image of the cosmic ray source in the following. When the cosmic ray source is fixed 
on the celestial sphere the image is the variation recorded during a sideral day. If, 
on the other hand, the cosmic ray source is related to phenomena within the solar 
system and gives an intensity variation with solar time the image is the variation 
recorded during a solar day. In the latter case the source moves on the celestial 
sphere. The aim of the present paper is to show how the image is influenced by the 
geomagnetic deflection of the primary particles and by the different processes 


Fig. 1. The extension of a cosmic ray source on the celestial sphere is the region enclosed by the 

curve on the left of the figure. The particle intensity along the dashed strip of the source is re- 

corded by the telescope when the earth rotates. @ 00 is the measuring solid angle of the telescope 
outside the earth’s magnetic field. 
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occurring in the atmosphere. The optics of cosmic ray telescopes thus includes a 
study of these two main parts as well as the influence of the geometrical dimensions 
of the telescope itself. 

In order to find the image of a cosmic ray source we start by finding the image of 
a source element subtending an infinitely small solid angle. Using the results of this 
study it is then possible to deduce the image of any extended cosmic ray source. 

The problem is divided into two main parts: a study of the deflection of the primary 
particles in the geomagnetic field and the influence of the atmosphere upon the 
particle flow. Variations in atmospheric conditions which introduce variations in 
cosmic ray intensity are, however, not considered in this paper. 


1. GEOMAGNETIC DEFLECTION OF PRIMARY COSMIC RAY 
PARTICLES 


Our knowledge in this field is extensive (although not yet sufficient) thanks to the 
mathematical calculations by Stérmer [1], Lemaitre and Vallarta [2], Firror [3], 
Jory [4], List, Schliiter and Katterbach [5], List [6] and many others [cf. 7, 8]. 
Scale model experiments by Malmfors [9] and Brunberg and Dattner [10] have also 
given important information about the geomagnetic deflection. 

In the calculations and experiments mentioned above, only the dipole term of 
the earth’s magnetic field has been considered. Terms of higher order may necessitate 
small corrections, principally in the lower momentum range, but little is known 
about these corrections at present. 


1.1. Discussions of the primary particle flow 
Consider a very small source of cosmic ray particles (a source element) situated 
far away from the earth, i.e. outside the region where the geomagnetic field has any 
deflecting influence on the particles. 
Let the source element subtend the solid angle dQ and dA be a surface element in 


the cross-section of the beam from the source element. 
The differential momentum spectrum is given by: 


I =I(p) particles/em?-sec:sterad - momentum interval. 


In the six-dimensional phase space the volume element occupied by beam particles 
with momenta between p and p + dp entering during the time dt is: 


21z 
dV,= (? aaa) (p2 dp dQ). (1) 


m denotes the relativistic mass of a particle. 
The number of particles within the element d V, is: 
dN =I(p)dpdA dQ dt. (2) 


Consequently the density of particles in phase space is: 
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_aN I(p)dp 


aed is p 


Liouville’s theorem of statistical mechanics shows that the density of particles 
within the element dV, is constant when the element moves about in phase space. 
This has been pointed out by Lemaitre and Vallarta [11] and further treated by Swann 
[12]. Thus particles originating from a surface element dA of the primary beam 
form a tube of flow in phase space with a particle density given by eq. 3. Compare 
also Janossy [13, p. 266]. 

As dA and dQ are infinitely small the tube of flow is infinitely thin. The flow 
pattern is governed by the geomagnetic deflection. 


D m. (3) 


1.2. The element in phase space representing a telescope 


We consider a telescope on the earth with infinitely small solid angle dw around 
pr and surface da perpendicular to p;, where the vector p, represents the measuring 
direction of the telescope, see fig. 2. 

The six-dimensional element dV; in phase space contains the number of particles 
with momenta between p, and p; +d, which are recorded by the telescope during 
the time dt. In the following d V7 is called the telescope element in phase space and is 
given by: 


dV >=" dtdaprdprdo. (4) 


Fig. 2. The projection of the solid angle w of a telescope located on the earth. The direction of 
the axis of the telescope is defined by pz. The projection is carried out by means of particle orbits 
in the geomagnetic field and the projected solid angle woo is visualized on the sphere, 
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The position of dV; in phase space changes with the direction of p; and is conse- 
quently dependent both on time (as the earth rotates) and the measuring zenith and 
aximuth directions of the telescope. When the telescope element is completely inside 
the particle flow in phase space, as regards both its momentum and three dimensional 
space parts, the number of particles within the element is: 


dn = DdVy, (5) 
where D is given by eq. 3. 
As the particle energy is constant in the magnetic field we have p = pz. From this 
and equations (3), (4) and (5) the number of particles recorded by the telescope per 
unit time is: 


dy =dn/dt=I(p)dpdadw particles/sec.(?) (6) 


From eq. 6 we conclude that a telescope measures the same intensity in the geomag- 
netic field as another telescope would measure outside the geomagnetic field, provided 
that the particle flow covers the surfaces and solid angles of the telescopes and that 
the product of the solid angle and surface of the two telescopes are equal. This means 
that: 

dadw =da,, do, (7) 


where da,, and dw,, are the surface and solid angle respectively of a telescope outside 
the geomagnetic field. 

The condition that the telescope element dV, must lie completely inside the flow 
in phase space for eq. 6 to be valid is interpreted in the following way: 

The momentum and three-dimensional space parts of a primary particle element 
dV, are respectively (compare eq. 1): 


AV pm =p? dQ dp, (8a) 
dV,,= dtd A. (8b) 
Ds m 


The motion up to the point in phase space where the telescope element d V 7 is situated 
transforms the two space parts into: 


(p? dQ dp),, (9a) 
p 
(2 aa4),. (9b) 


The transformation is governed by the particle motion in the geomagnetic field, as 
will be discussed in section 1.3. The telescope records those particles which are within 
those parts of the transformed elements of equations 9 which coincide with the corre- 
sponding parts of dV,. Full intensity is measured when the transformed elements 
completely coincide with or extend beyond the corresponding parts of dV: 


(p? dQ dp), > py? dw apr, (10a) 


Pp Pr 
—dtdA})r>— 
(* ): a dtda, (10 b) 
(1) The number 3 denotes the order of the differential. 
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Or as Pr = Pp: dQ, >dw, (1la) 
dA, >da. (11b) 

The condition of coincidence naturally requires: 


1) That the whole range of directions of particle motion at the telescope com- 
pletely covers or extends beyond the solid angle of the telescope. 

2) That the transformed surface element of the cross-section of the beam com- 
pletely covers or extends beyond the surface of the telescope. 


As will be shown later it is more convenient from the practical point of view to go 
the other way and transform the telescope element from the earth to infinity. A 
discussion similar to the one above then gives instead of equations lla and 11b: 


dQ > daw, (12a) 
dA >daq. (12b) 


The condition of coincidence is then that the whole range of directions of particle 
motion in the primary beam completely covers or extends beyond the transformed 
solid angle dw,,. Eq. 12b is always fulfilled as the cross-section of the primary 
beam was assumed to be very large. 

We thus summarize that a telescope with infinitely small solid angle and surface 
records full particles intensity from the momentum interval dp in a beam of large 
cross-section subtending the solid angle dQ2 when the condition of coincidence is 
fulfilled together with the equations 7, 12a and 12b. The counting rate is derived from: 


d'y = I(p)dpdadw =I (p)dpda,dw,..=I(p)dp da (dw/dw~)do., (18a) 
and full counting rate is: 
d*y = I(p)dp da (dw/dw..) dQ particles/sec. (13b) 


The telescope measures either full intensity or nothing, when considering particles 
in a given momentum interval. Full intensity means that the telescope element in 
phase space is within a tube of flow, which happens when the direction of dw,, 
coincides with the direction of dQ. Zero intensity is measured when the telescope 
element does not lie within any tube of flow. 

We call the ratio between the infinitely small solid angle of the telescope and its 
projected solid angle the magnification M: 


M =dw/do.. (14) 


For infinitely high momenta it is evident that M—=1 but for lower momenta we 
may have M > 1 as well as M< 1. For certain momenta and geomagnetic latitudes 
M approaches infinitely high values which is connected with a focusing effect in 
the geomagnetic field. This has been discussed by Astrém [14]. At momenta near 
geomagnetic cut-off M approaches zero.(?) 


(4) This was clearly observed when determining the geomagnetic deflection experimentally, 
Brunberg and Dattner [10]. The luminous spot on the fluorescent screen, which for high momenta 
was approximately 2 cm*, increased enormously near the cut-off momenta. The surface of the 
spot could be used for an accurate determination of M. 
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We have, from equations 12b, 13a and 14: 
dA =da (dw/dw,.) =da M, (15) 


ie. the surface element d.A which completely coincides with da, depends on M and, 
consequently, so does the counting rate when the primary beam subtends an infi- 
nitely small solid angle but has infinite extension.(?) 

In the case when particles are moving isotropically in space around the earth, a 
telescope always measures the same intensity because the particle flow covers the 
projected surface and solid angle of a telescope for any allowed transformation through 
the geomagnetic field. This was first pointed out by Lemaitre and Vallarta [11]. 


1.3. The transformation of solid angle and surface of a telescope 


The transformations or projections are carried out along the particle orbits in the 
geomagnetic field. The direction of orbits far away from the earth are called asymp- 
totic directions. They are represented by vectors emanating from the centre of a globe 
and the direction of a vector is defined by the geographic latitude and longitude angles 
® and ‘’. The longitude angle is reckoned east from the meridian plane through the 
point where the particle orbit hits the earth. On the globe surface the end points of 
the vectors give the asymptotic directions for different particle momenta, and measur- 
ing directions on the earth. End points corresponding to different particle momenta 
but to the same measuring direction are joined, giving a curve on the globe. The 
curves for different measuring directions are used to find the projected solid angle 
of a given telescope as described by Brunberg [15]. 


Fig. 3, The sphere of asymptotic directions. p 
denotes the angular position of the source, 
from which particles of momentum p arrive, 
within the small solid angle AQ. The area 
on the sphere enclosed by the dashed line 
represents the projected solid angle woo of a 
S telescope on the earth. 


(*) This is in accordance with the general optical law of Helmholtz which, applied to particle 
motion in electric and magnetic fields, states that: 


ve cos 8) d Fy dQ, = vf cos #, dF, dQ,, 


where v = the particle velocity, d F = a surface element of the particle beam, dQ) = the solid 
, ’ = angle 
subtended by the beam, and # = the angle between the normal of the surface, d F, and the eat 


tion of the beam. Subscripts refer to different positions along the b 
lagen der Elektronenoptik, p. 623), ; : urbimdenentbips se 


In our case we have v, = v and cos # = 1. Consequently: 
dF, /d F, = dQ,/dQ, = M. 
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Fig. 2 shows the projection of a finite solid angle of a telescope for a given par- 
ticle momentum. In the following such a globe will be referred to as a sphere of 
asymptotic directions. The coincidence of the primary particle flow and the projec- 
ted finite solid angle of a telescope is easily visualized on the sphere, as shown 
in fig. 3, and is further treated in section 1.4. 


1.4. The properties of a telescope of finite dimensions 


In the preceding sections the recording properties of a telescope with infinitely 
small solid angle and surface have been discussed. Consider now a cylindrical tele- 
scope bounded by two circular surfaces, a, separated by a distance, d, as in fig. 4. 


1.4.1. The general case 


The counting rate of the telescope is (from eq. 6): 


dyv=I(p)dp{ | dadw. (16) 


Eq. 16 means an integration in phase space over all telescope elements dV, which 
are completely within the particle flow and thus contribute to the counting rate. Let us 
first consider the integration over the solid angle in eq. 16 and find the intensity 
recorded by a surface element da and the whole solid angle w, of the telescope belong- 
ing to da. 

Using equations 13a and 14 we have the following expression for the number of 
particles recorded by a surface element da of the given telescope: 


dy=I(p)dpda i cosuMda,, (17) 


(®o) 60 


where (@ )<. is the projection of the solid angle w, on the sphere of asymptotic direc- 
tions. (da cos u) is the projected area of da on a plane perpendicular to the direction 


Vertical 


2 7 Parallel 
a TS) ee 


Fig. 4. A cylindrical telescope bounded by two 

circular surfaces separated by a distance d. 

The shaded area is the recording part of the 

surface for particles arriving at an angle wu to 
the axis of the telescope. 
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of particles arriving at the telescope. Contributious to the integral come only from 
the part of the solid angle which contains sources; consequently: 


{i cos u Mdw..= {cos u M dQ, 


(®o) 66 


where ¢ is the extension of the source within (wo)... Thus the integration over the 
solid angle of the telescope is replaced by an integration over the solid angle part 
of the source which is within the projected solid angle of the telescope on the sphere 
of asymptotic directions. The counting rate due to a source element d{ within the 
projected solid angle is: 


dy = I(p)dpdacosu M dQ. (18a) 


As we shall discuss the counting rate of a telescope due to a source subtending a 
small solid angle we introduce in the following, for the sake of simplicity, a source 
element subtending the small but finite solid angle AQ. As M and u are assumed to 
be constant over the source, integration of eq. 18a gives: 


d*y = I(p)dpdacosu M AQ. (18b) 


Fig. 5a shows the situation. Three surface elements in a telescope are drawn as 
an example with their corresponding solid angles w,, w, and w,. The corresponding 
solid angles on the sphere of asymptotic directions are shown in fig. 5b. The dashed 
curve is the limit of projected solid angles belonging to surface elements on the 
fringe of the lower telescope surface. The total solid angle of the telescope is defined 
by the apex angle 2 Uma, with Umax =are tg(2r/d) and equals 27(1 — cos Umax), see 
fig. 4. The position of AQ on the sphere indicates the angular position of the source 
element and determines the angles of incidence at the telescope. u is the angle between 


Fig. 5. (a) Three surface elements with corresponding solid ang] i 
(a) 1 gles are chosen as an example in 
the cylindrical telescope. (6) The solid angles projected on the sphere of asymptotic dientsaran 
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the central telescope direction and the direction of arriving particles. According to 
equation 18b the surface elements da, and da, contribute to the counting rate by: 


d*y, = I(p)dp da, cosu M AQ, (19a) 
d*y, = I(p)dp da, cos u M AQ, (19b) 


whereas there is no contribution from dag. 
The surface integral yields: 


dy =I(p)dpacosu G(u) M AQ. (20) 


G(u) is a function determined by the geometry of the telescope and takes into 
account the fact that only a part of the surface records particles arriving from the 
direction defined by the angle wu. This part is shaded in fig. 5a where the surface 
element da; is outside the recording region. As mentioned above, wu is a function of 
the positional angles ® and of the source element on the sphere of asymptotic 
directions. The function cos u G@(u) is calculated in appendix I for cylindrical tele- 
scopes of different apertures. As the magnification M depends on the orbits in the 
magnetic field it is a function both of the particle momentum and the angular posi- 
tion of the source element. 

The sphere of asymptotic directions is used when calculating the counting rate of 
a telescope when the earth rotates relative to the source element. The total solid 
angle of the telescope is projected on the sphere using particle orbits corresponding 
to the momentum » as in fig. 3. The direction to the source element is given by the 
vector p and the source subtends the solid angle AQ. When the earth rotates AQ 
moves along a parallel circle on the sphere. In the angular interval AY’ the counting 
rate of the telescope is determined by eq. 20. In this way it is possible to calculate 
the counting rate due to different parts of the source spectrum for any angular posi- 
tion (®,‘f’) of the source element. The total counting rate of the telescope is obtained 
by adding the contributions from all parts of the spectrum, 


1.4.2. The narrow angle telescope 


For small values of «, i.e. for telescopes with the separation between the surfaces 
large compared to the radius of a surface, we have: 


cos uw 1, M ro/ox, wo wa/d. 
The counting rate is from eq. 20: 


dv~1(p)dp AQ a Gu), (21) 


1.4.3. The ionization chamber 


In a spherical ionization chamber the magnitude of the measuring surface as well 
as the solid angle is the same in all directions. The chamber is said to be omnidirec- 
tional. The measuring surface (cross-section through the centre of the chamber) 
perpendicular to the direction of arriving particles shall of course be regarded as a 
fictious one as the particles ionize along their whole path within the chamber. Particles 
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passing through the centre of the surface consequently contribute more to the re- 
corded intensity than other particles. But this variation of sensitivity over a cross- 
section of the chamber does not complicate the calculations as it is the same in all 
measuring directions. The counting rate according to eq. 20 with G(u) =1 is: 


dy =I(p)dp AQ Mag, (22) 


where a, is the effective measuring surface of the chamber. 

When calculating the counting rate of the chamber, its solid angle, 27 (only the 
upper half sphere) is projected on the sphere of asymptotic directions. As long as the 
angular position of the source element is within the projected solid angle the counting 
rate is given by eq. 22. 


1.4.4. The three-dimensional polar diagram of a telescop2 on the sphere of 
asymptotic directions 


The counting rate of a telescope is calculated according to section 1.4.1 for different 
angular positions of the source element represented by the solid angle AQ on the 
sphere of asymptotic directions and assuming a certain momentum spectrum. Each 
position of the element gives a certain counting rate in the telescope for a given 
momentum interval—the sensitivity of the telescope varies with the position of AQ. 
At every position of AQ we mark out a distance in the radial direction proportional 
to the counting rate of the telescope with AQ in this position. This gives a curved 
surface above the globe surface, the height of which varies continuously over the 
globe. This ‘‘mountain”’ is the three-dimensional polar diagram of the telescope in 
the given momentum interval. The total three-dimensional polar diagram is suitably 
found by dividing the whole spectrum into a number of finite momentum intervals 
Ap and calculating the curved surface given by each interval. The geomagnetic de- 
flection of particles within an interval is then represented by the mean value of p in 
that interval. Superposition of the contributions from the whole spectrum constitutes 
the total three-dimensional polar diagram of the telescope on the sphere of asymptotic 
directions. 

In the following, a section along a parallel circle of the curved surface will be used. 
This represents the counting rate of the telescope as a function of ’, i.e. the longitudi- 
nal position of the source reckoned to the east from the meridian plane through the 
place of observation on the earth. This section is important in the following discus- 


sions and is referred to as the polar diagram of the telescope along the parallel circle 
in question. 


1.4.5. A few examples of calculating polar diagrams when the atmospheric influence 
is neglected 


We shall now discuss the calculation of polar diagrams of an ionization chamber 
and of cylindrical telescopes with different apertures but with the same surfaces. 
The telescopes point in the zenith direction and their geometrical configuration is 
shown in fig. 4 where now £, = 0. 

The angular position of the source is along a parallel circle on the geographic 
latitude ® and the beam from the source subtends the small solid angle AQ. The 
source spectrum used is the same as that of the isotropic background which seems 
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I(p) particles/s cm* sterad momentum unit. 


Fig. 6. (a) The primary differential mo- 

mentum spectrum schematically drawn. 

(6) The positions of the projected solid 
angle for different primary momenta. 


Direction of the beam 
from the source. 


reasonable in this first approach. The differential spectrum adopted is proportional 
to p *, ef. [46]. 

The total solid angle of the telescope which corresponds to the apex angle 2 wax 
is projected on the sphere of asymptotic directions. This is done for different momenta, 
each of them representing a certain momentum interval dp. The projected total solid 
angle on the sphere contains all possible directions of incidence at the telescope for 
momenta between p and p+dp. When the angular position of the source element 
is within the projected solid angle the counting rate due particles in the momentum 
interval dp is given by eq. 20. 

For a certain angular position of the source element the total counting rate is: 


Pmax 
y=aAQ | M cos u G(u)I(p) dp. (23) 


Pmin 


Poin ANd Pmax are the minimum and maximum momenta for which the corresponding 
_ projected solid angles contain the source element on the sphere of asymptotic direc- 
tions. This is shown schematically in fig. 6. The counting rate is computed for a 
number of source positions along the parallel circle ® = const. and the counting rate 
thus obtained as a function of the longitudinal angle VY is the polar diagram of the 
telescope along this parallel circle. 

In the following the spectrum has been divided into a number of finite intervals 
Ap with each interval centred at a momentum p;. The number of particles within 
that interval, I(p,;) Ap, are taken as having the same deflection represented by the 
momentum p, and their asymptotic directions are within the corresponding solid 
angle on the sphere of asymptotic directions, fig. 7a and 6. 

For each position of AQ on the sphere the counting rate due to particles in the 
momentum interval between p; and (p; + Ap) is: 


Ay, =a AQ I(p,;) Ap M, cos u; G(x). (24) 
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I (p) particles/s cm? sterad momentum unit. 


c. Fig. 7. The principle of calculating a pola- 
diagram along a parallel circle. (a) The pri- 
mary differentialmomentum spectrum. (6) The 
projected solid angle for different momenta. 
(c) The polar diagram along a parallel circle. 


particles /sec. 


The magnification M and the angle of incidence («,f) at the telescope are found from 
the sphere. As regards the angles, only the zenith angle 6 matters in this exampel 
because the central direction of the telescope points in the zenith direction, from 
which it follows that «=f. The counting rate is thus written: 


Ay, =a AQ I(p,) Ap M, cos B; G(8)). (25) 


cos 6 G(B) has been calculated in appendix I. 
The counting rate of an ionization chamber is (compare eq. 22): 


Ay, =a, AQ I(p;) Ap M;. (26) 


The counting rate is computed in this way for a number of angular positions ¥’ of 
the source within the total solid angle and plotted as a function of V, see fig. 7), c. 
The polar diagram thus obtained (for instance P(p,) in fig. 7c) is approximately that 
corresponding to particles in the momentum interval between (p,—}Ap) and 
(p, +} Ap). 

The approximation is naturally better the smaller Ap, but as the accuracy of the 
projected solid angles amounts to 4°-5° in position, sometimes more, this puts a 
lower limit to Ap. The limit value could, however, be increased still more as, in 
reality, the measured counting rate is the averge taken over a considerable time, 
usually two hours, which corresponds to a change of 30° in the position of AQ on 
the sphere. 

The procedure described above is then repeated for other momentum intervals and 
the total polar diagram is obtained by adding up the contributions from the whole 


spectrum, fig. 7c. Now the relative displacement of the polar diagrams corresponding 
to different momentum intervals puts an upper limit to Ap. 
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Fig. 8. An example of calculating polar diagrams along ® = 0 in the momentum range around 


50 GeV/c at Huancaoy. 

(a) The angle of incidence f at the telescope for different angular positions of the source along 
the equator. (b) The magnification M. (c) Polar diagrams of zenith pointing telescopes with 
different apertures, wmax. 


Fig. 8 illustrates the method of calculating polar diagrams in the momentum 
interval centred at 50 GeV/c at the Huancayo station (geomagnetic latitude = 0.6° S) 
with the source element in the equatorial plane (0 =0). In fig. 8a the zenith angle 
of incidence, f, at the observation place is plotted. 6 is found from the Huancayo 
sphere of asymptotic directions by the different positions ‘’ of the source, cf. [15]. 
Fig. 8b shows the corresponding magnification M as a function of ¥. 

The magnification has been approximately determined by projecting the solid 
angle of a circular telescope with wma, = 8° on the sphere along the equator for dif- 
ferent directions of incidence, corresponding to the positions of the source. The 
ratio Ws:/(Ws-).o is then taken as representing the magnification M. The accuracy of 
M when determined in this way is of the order of 0.2—0.3 in the chosen example but 
for lower momenta the error is considerable and it is finally impossible to determine VM. 
This is due to the fact that the projected solid angles are deformed and their size 
increased enormously when approaching cut-off momenta. As for each zenith angle 
only four asymptotic directions are known (N, E, 8 and W directions) it is thus 
impossible at present to find the projected solid angles near cut-off momenta on the 


sphere. 
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Fig. 9. (a) The contributions from different momentum invervals to the polar diagram along 
® = 0 of an ionization chamber in Huancayo. (6b) Polar diagrams of different telescopes when the 
primary spectrum in the range 20-100 GeV/c has been used. 

Dot-dashed lines represent the probable shape when the contributions from momenta below 
20 GeV/c have been included (rough extrapolation). The differential momentum spectrum adop- 
ted is proportional to p~?. 


Fig. 8c shows the polar diagrams obtained with different telescopes using equa- 
tions 25 and 26 where I(p) Ap AQ is for convenience put equal to one. 

Fig. 9a shows the polar diagram along ® = 0 of an ionization chamber when the 
contributions from the primary spectrum in the range 20-100 GeV/c have been 
added with regard to the proper weights given by the momentum spectrum. The 
different areas of the polar diagram represent the contributions from different mo- 
mentum intervals. Dashed lines mean that the values are uncertain which is mainly 
due to the errors in M as pointed out above. It is seen that the low energy part of 
the spectrum plays an important role for the shape of the polar diagram. It is conse- 
quently necessary to calculate more particle orbits in this region. 

In fig. 96, finally, the polar diagrams of cosmic ray telescopes with different aper- 
tures are drawn. These diagrams are of great importance for finding the image of 
a cosmic ray source with the telescope situated above the atmosphere at the Huan- 
cayo station. They are discussed further in chapter 3 where the influence of the 
atmosphere has also been taken into account. 


2. THE INFLUENCE OF THE ATMOSPHERE 


A cosmic ray telescope at sea level or deep down in the atmosphere records almost 
only secondary particles produced by nuclear interactions in the atmosphere. The 
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atmospheric influence upon the polar diagram of a telescope is twofold. First the 
contributions from different parts of the primary momentum spectrum to the polar 
diagram are not the same as when measuring at the top of the atmosphere, because 
the production rate and energies of the secondaries are functions of the primary 
particle energy. Secondly several processes occur in the atmosphere which have a 
spreading effect upon the secondaries so that a parallel beam of primary particles 
is gradually converted into a more or less divergent beam of secondaries which are 
finally recorded by the telescope. This makes the polar diagram wider. 


2.1. The spread in the atmosphere 


In the beginning of cosmic ray research it was assumed as a first approximation 
that the direction of particle motion is conserved in the atmosphere. In 1951 Vernov 
et al. [16] deduced experimentally an expression for the angular distribution of sec- 
ondaries. They found a mean angle of deviation of about 30° between the directions 
of the secondary and primary particles, i.e. about half the number of the ionizing 
secondaries are found within a cone with the apex angle 60°. Sandstrém [17] has 
pointed out that the spreading effect of the atmosphere should be taken into account 
when constructing cosmic ray telescopes. 

In this chapter the spread of a parallel beam of primary particles is estimated. The 
different processes occurring in the atmosphere are not completely known; for example 
the distribution of mesons produced in nucleon-nucleon collisions at different primary 
energies. In such cases we must make certain assumptions in order to be able to carry 
the calculations to the final stage and be able to survey the field. 

Only the meson component of the secondary flow is considered; that is, we neglect 
the electron, photon and nucleonic components of the particle flow. Furthermore 
we assume that the recording apparatus is shielded by 10 cm of lead thus giving re- 
sponse only to w-mesons of energies higher than 0.1 GeV. 

We also assume that a primary particle, after having performed one nuclear interac- 
tion with an air nucleon, does not give any further contribution to the ~-meson 
flow or to the counting rate of a telescope at sea level. 

We shall consider three parts of the problem: 

(1) The rate of inelastic nucleon-nucleon collisions throughout the whole atmosphere. 

(2) The intensity and energy distribution of mesons created in a nucleon-nucleon 
collision. 

(3) Processes influencing the secondaries on their way towards the recording ap- 
paratus. 


2.1.1. The rate of inelastic nucleon-nucleon collisions in the atmosphere 


We consider a nearly parallel beam of protons with the energy W,, entering the 
earth’s atmosphere from a direction making the angle f’ to the zenith direction, as 
in fig. 10. The probability of interaction between a primary proton and an air nucleon 
per unit path length depends on the mass density along the path. If 7) particles/em*- 
sec:sterad. enter the atmosphere in the beam, the number of particles at the height 
z which have not undergone nuclear interaction is: 


Nz pr = Ng EXP (— Mzp/z)- (27) 
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Fig. 10. A telescope at sea level with the infinitely small solid angle dw records mesons arriving 
in the direction «, 8. The measons are produced at the height z by primary particles coming from 
the direction «’, B’. 


mzg is the penetrated mass of air in g/cm? from the top of the atmosphere along 
a path which makes the angle f’ with the vertical. L is the mass corresponding to 
the mean free path of nuclear interaction for primary particles in air. 

Next we want to know the number of particles colliding inelastically, i.e. undergo- 
ing nuclear interactions between the heights z and z + dz in the atmosphere which 
is static and isothermal and thus described by the variation of density: 

g=Q0e", (28) 
where | = mo/0 = the scale height, assumed const.; my =the mass per cm? at sea 
level; 09 = air density at sea level. 


The penetrated mass from the top of the atmosphere to the height z is: 


m “7 
0 e zit 


Map = oF (29) 


The number of inelastic collisions between the heights z and z + dz is from equa- 
tions 27 and 29 by derivation: 


z 


dn AP asi asa EARLE do d ticl : d 
hl ihad FY B P 1 Leos Bp’ p 7) |¢2 partic es/cm?- sec: sterad. (30) 
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Fig. 11. The number of inelastic collisions in the atmosphere as a function of the height for 
different angles of incidence to the vertical. 


dn-g is consequently the change in the number of primary particles of the beam 
which have not undergone nuclear interaction, in a layer dz of the atmosphere.(') 
The number of inelastic collisions per unit volume is (from eq. 30): 


Bes, mo ex oye wes a ex pare articles 3 terad Di 
dz ° Li cos p’ Puleitg L cos p’ P ] particles/em*-sec-sterad. (31) 


The number of inelastic collisions per unit volume according to eq. 31 has been 
plotted in fig. 11 as a function of the altitude for different values of #’ and with 


(1) The atmospheric depth in g/em* taken from the top of the atmosphere is often used in 
expressions of this type instead of the height z above sea level. As we are interested in knowing 
where the production rate of secondaries is a maximum we are forced to introduce the mass 
distribution of the atmosphere. Furthermore, the path length of secondaries enters into the caleu- 
lations, which is another reason for using the height z in the present equations. 
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Fig. 12. The mass distribution in the atmosphere. 


Mm, = 1030 g/cm? and 1 = 8000 m. L probably depends little on the primary energy 
and experiments give values between 100-150 g/cm? (compare Rossi [18, p. 442], 
and Sands [19]). The value 120 g/cm? has been used in the present calculations. 

We define the production rate of secondaries as the product of the number of nuclear 
interactions and the average number of mesons produced per interaction. 

The maximum production rate for primary particles arriving in the zenith direc- 
tion occurs at a height of 17 km, as seen from fig. 11. This corresponds to a pressure 
of 120 g/cm? with the assumed atmosphere or to 90 g/cm? using the ‘““NACA-atmos- 
phere” (National Advisory Committee for Aeronautics), fig. 12. 

The number of inelastic collisions per unit volume, integrated over all primary 


directions, assuming isotropic distribution of particles over the hemisphere and above 
the atmosphere, is: 


22x x/2 
dn _ Ny Mo _2 My z RY ee 
dz | x aes I Ae ‘| Saas ee 
LA 
particles/em-sec (32) 
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Fig. 13. Curve A: The total number of inelastic collisions per unit volume in the atmosphere as 

a function of the height. Curve B: The number of primary particles from the zenith direction 

which have not collided inelastically. Curve C: The total number of primary particles which have 
not collided inelastically (particles/em? - sec). 


where a’ is the azimuth angle of the primary beam, as in fig. 10. The integration has 
been carried out in appendix II and the total number of inelastic collisions per unit 
volume plotted in fig. 13. The maximum total production rate according to these 
calculations occurs at a height of around 24 km, which corresponds to a pressure of 
30 g/cm? in the NACA-atmosphere. 

Measurements of the ionization in the atmosphere with balloons and rockets using 
unshielded ionization chambers show a maximum at a pressure between 40-110 g/cm?, 
depending on the geomagnetic latitude 2. The lower value is at A = 60° N and the 
high one near the equator according to Bowen, Millikan and Neher [20], Gangnes 
et al. [21]. 

The vertical intensity measured by Vidale [22] has a maximum at 75 g/cm? 
(2 = 55°) and 135 g/cm? (A = 28°) for the total radiation, whereas the corresponding 
figures for the hard component are approximately 48 and 120 g/cm?. The maximum 
ionization is thus measured lower down in the atmosphere compared to the level of 
the maximum production rate (30 g/em*) found in our calculations. 

This is however what would be expected since an instrument records not only sec- 
ondaries produced near it but also secondaries coming from higher up in the atmos- 
phere as well as the primary particle intensity itself. If all secondaries could reach sea 
level the intensity would obviously increase with decreasing altitude and highest 
intensity be measured at sea level. If, on the other hand, the number of secondaries 
produced at a certain altitude were to decrease with a very short mean free path, 
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maximum intensity would be measured near the level where the production rate is 
maximum. In reality there is a decrease in intensity due to absorption and decay which 
shifts the maximum towards lower altitudes.(') 

The observed latitude dependence of the height of maximal ionization is probably 
due to the lower mean energy of secondaries which are produced by the low energy 
part of the primary spectrum. As these secondaries decay more rapidly and are ab- 
sorbed after a shorter path length the maximum ionization should be measured 
nearer to the level of maximum production, which according to our calculations is 
at 30 g/cm?. The measured shift of the level of maximum ionization from 110-130 
g/cm? at lower latitudes up to 40-80 g/cm? at high latitudes, where there are more 
primaries of low energy, is in qualitative agreement with the hypotheses of a height 
difference between the levels of maximum production rate and maximum ionization. 

When discussing the production rate in the atmosphere it should be pointed out 
that the zenith angle dependence of the intensity at different altitudes is probably 
explained by the number of inelastic collisions per unit volume given by eq. 31 or 
properly speaking by the production rate. Different authors have found that the 
intensity increases with zenith angle above a certain altitude (Winckler and Stroud 


(*) For example the decay of u-mesons with energy Wy, is given by: 
Nu = Now exp(—s/cT), 


where s = the path length; t = t’ (Wy/myc*); and 7’ is the decay time of a p-meson at rest. With 
Wu = 1 GeV the mean free path of decay is 6000 m. 
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Fig. 15. Zenithal and azimuthal distribution of cosmic ray intensity at the equator (Vernov et al.). 
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[23], Vernov et al. [24]). It has been suggested by Singer [25] that this is due to an 
albedo effect. 

Winckler and Stroud, however, have compared the measurements of the zenithal 
particle distribution with their own calculations of the expected distribution, assum- 
ing different absorption coefficients for the primary and secondary particles. They 
have investigated the influence of different angular distributions of secondaries 
and have found that the closest agreement is reached by assuming a sharply colli- 
mated forward type of distribution. 

From eq. 31 the number of inelastic collisions per unit volume in the atmosphere 
has been calculated as a function of the zenith angle /’ with the height to the produc- 
tion level as a parameter. This has been plotted in fig. 14. Assuming as a first approxi- 
mation that all the mesons produced in nucleon-nucleon collisions move in the same 
direction as the arriving protons, fig. 14 gives the production rate of secondaries in 
different directions in the atmosphere. We notice that up to about 15000 m the produc- 
tion rate is a decreasing function of /’, at 17000 m it is constant for Bp’ <45° and 
above 20000 m it increases with the zenith angle. 

Allowing for a shift towards lower altitudes according to what was discussed above 
these curves are in good agreement with the measurements of Vernov et al. and 
Winckler and Stroud (compare figs. 14, 15 and 16). We thus conclude that eq. 32 and 
figs. 12, 13 and 14 are representative of the production rate in the atmosphere. 

Finally, assuming isotropic particle distribution over the hemisphere above the 
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Fig. 16. The zenithal distribution of cosmic ray intensity measured at 56° geomagnetic latitude 
(Winckler and Stroud). 


atmosphere, the zenith angle dependence of the intensity at sea level could be derived 
by superimposing the curves of fig. 14, with the decay and absorption of “-mesons 
taken into consideration. This means that the contribution from different levels 
decreases with increasing height. A calculation along these lines is carried out in 
section 2.3. 


2.1.2, The intensity and energy distribution of mesons produced in nucleon-nucleon 
collisions 


The interaction between a primary cosmic ray particle and an air nucleon results 
in the production of secondary particles, as for example other nucleons, electrons and 
different kinds of mesons. For the primary energies which are most interesting in 
our problem the production of z-mesons is important. We want to know the intensity 
and energy distribution of the produced z-mesons as function of the angle between 
the direction of the arriving primary particle and the direction of the emitted mesons 
in a coordinate system at rest (Z-system). 

This is done firstly by finding the intensity and energy distribution of the mesons 
in the centre of mass system of the two colliding nucleons (C-system). Secondly, these 
distributions are transformed to the L-system. 

The intensity distribution of the produced mesons in the C-system is probably 
anisotropic at high primary energies, as pointed out by Schein e¢ al. [26], but in the 
range around 10 GeV the distribution is more likely to be isotropic. This is assumed 
in the following calculations in the absence of more precise information [27, 28]. 

The energy distribution of z-mesons in the C-system has been studied by Camerini 
et al. [29]. Under the assumption that all the mesons are produced in a single nucleon- 
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Fig. 17. (a) The angular intensity distribution of a-mesons in the L-system. Each curve corre- 

sponds to z-mesons produced isotropically in the C-system all with the same energy. In the 

energy case corresponding to m = 1.10 there are no mesons in the L-system for § > 33°. This is 

explained by Bradt, Kaplan and Peters [30]. (b) The angular energy distribution of 2-mesons 
in the L-system. 


nucleon collision these authors find a differential spectrum which is approximately 
represented by: 
dnz V(y,)2—1 
RA ies (33) 
dy, (Ya) 


where y, is the energy of a z-meson in C-system in terms of its rest energy. 

The transformation of the intensity and energy distributions to the L-system has 
been treated in great detail by Bradt, Kaplan and Peters [30]. The complete formulae 
for the angular intensity distribution J (0, W,) and the z-meson energy expressed in 
units of the rest energy y,(0,W,) in the L-system are collected in appendix III. 

The shape of the distribution function /(0,W,) in the L-system is determined 
mainly by the primary proton energy W,, but also by m, i.e. the ratio between the 
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proton and meson velocities in the C-system. The value of meson velocity v, in the 
C-system is dependent on the inelasticity of the collision (i.e. the ratio between the 
energy of the produced mesons and the total energy available) and the total multi- 
plicity which includes both charged and uncharged z-mesons. 

The following calculations have been carried out for 20 GeV/c protons as an 
example, because this momentum seems to be representative of the primary spectrum 
when measuring the position of the cosmic ray source, which is responsible for the 
daily variation of cosmic ray intensity (Brunberg and Dattner [31)). 

The Camerini-distribution (eq. 33) has been approximated by a step function by 
choosing three energy intervals in such a way that the same number of mesons are 
contained in each energy interval. The mean meson energy of each interval, together 
with assumed primary energy W,, yields the values of m = 1.10, 1.00 and 0.96 with 
the corresponding values of y, = 2.0, 3.35 and 10.0. For these values the distributions 
of intensity and kinetic energy of z-mesons in the L-system have been calculated 
using the formulae of appendix III. The results are shown in fig. 17. 


2.1.3. Different processes influencing the intensity and energy distributions 
of mesons 


On their way down through the atmosphere the charged u-mesons interact with air 
nuclei or decay into “-mesons and neutrinos. The number of uw-mesons then decreases 
gradually due to absorption and decay. Finally, Coulomb-scattering and deflection 
in the earth’s magnetic field have a spreading effect on the u-meson beam. 


2.1.3.1. The ty decay 


The mean life of a z-meson before it decays into a u-meson and a neutrino is about 
2.5 x 10-8 s. The mean free path is calculated from : 


Cc. 


My C2” 


ly 


(34) 


where t’ is the decay time in the system at rest. 

With W, =3 GeV we get 1, =160 m which is negligible in the calculations. 

The decay of the z-mesons introduces a small spreading effect, so that the ~-meson 
beam is slightly more divergent than the original z-meson beam. The magnitude of 
this effect is estimated by calculating the angular half-width 6, in the L-system of 
the u-meson intensity distribution. The angular half-width 0, is defined by a cone 
with the apex angle 20, within which half the number of the produced j-mesons is 
found in the Z-system. We use the modified eq. AIII.3 in appendix III: 


sin 9’ / Un\? 
t: = ——<———— - ie st = 
Bo cos 0 + v_/ vy, : ("") (Se 


As, in the C-system (the moving z-meson), any direction is equally probable for the 
created #-meson, the distribution in the C-system of a large number of events is 
isotropic. The half-width in the L-system is thus obtained by putting 6’ =42 in eq. 35: 


v, v,\2 
t ~%y/ lol 
Be tha A ("") (36) 
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Fig. 18. The angular energy distribution of mesons 

in the L-system for m = 1. Curve A: The total energy 

of mz-mesons. Curve B: The mean kinetic energy of 
fu-mesons due to the 7p decay. 


Using the formulae for particle decay to determine the momentum p, ~30 MeV 
(cf. Heisenberg [7], p. 561) and 3 GeV as the z-meson energy as an example, the half- 
angle amounts to 0.9° which is quite negligible compared to the half-angle of the 
intensity distribution of z-mesons which, with eq. AIII.19 in appendix III, is found 
to be 18° for 20 GeV/c primary protons. 

The energy distribution is, however, modified by the decay. z-mesons of energy 
W,,, emitted at an angle 9 to the direction of the primary beam, form a spectrum of 
y-mesons of energies ranging from a maximum to a minimum value and with a small 
spread in angle as discussed above. This narrow interval can be looked upon as a 
compressed strip at the angle @ containing all energies between the maximum and 
minimum values, as pointed out by Ascoli [32]. 

The w-meson energy is given there by: 


1 [Ma , Mul My C2 je : 
Wir=s |e me] mee Wt BE py, cond (37) 
or W.=0.789 Watcp, V2 —1xcos 6’, (38) 


where p, is the momentum of the “-meson in the C-system due to the decay. 
With W, =3.0 GeV eq. 38 is written: 


W,, = 2.36 + 0.63 x cos 6’ GeV. (39) 


As any energy in the interval between (Wy)max and (W,)min is equally probable we 
give each meson the mean value: 

<W > = 0.789 W,. (40) 

The angular distribution of kinetic energy of the y-mesons is calculated in this 


way and plotted in fig. 18. 
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2.1.3.2. The decay and absorption of u-mesons 


As shown in the previous section, the angular distribution of meson intensity 
remains practically unchanged by the 71> decay. 

The angular distribution of -mesons originating from a thin atmospheric layer dz 
at the height z is, according to eq. 31 in section 2.1.1.: 


Ny m Z m z : 
dn,z=I1(60, Wp) Ll cos ’ ee exp |- af ae B’ exp (- ‘)| dzdw 
particles/em?+sec-sterad., (41) 


where dw’ is the solid angle subtended by the primary cosmic ray beam and J (0, W,) 
the angular intensity distribution of produced mesons in a nucleon-nucleon collision. 

If the -mesons were to travel without losses or any other disturbing effects, a tele- 
scope with the surface da and an infinitely narrow solid angle dw perpendicular to 
the surface would record: 


dn, =dadwdn,z  particles/sec. (42) 


Due to decay, the particle density of the u-meson beam decreases and a telescope 
at sea level measures: 


Bn, =dadwdn,, exp (- particles/sec., (43) 


ak Ee 
Unt Cos B 


where t = t’ (W,,/m,c?) and v, = the velocity of the u-meson which is approximately 
equal to that of light. 

Introducing the numerical values: t’ = 2.15 x 10~* sec., m,c? =0.107 GeV, vz = 
3 x 108 m/s, eq. 43 is written: 


Bn, = dadw Any 2 exp Geary), 5 (44) 
i" a 


where C, = 60.3 x 10? if W,, is expressed in units of 10° volts. The change of W, along 
the path, due to atmospheric losses, modifies the distribution given by eq. 44 but 
this effect is found to be negligible. 

The “-mesons lose energy through ionization and excitation in the atmosphere. 
The magnitude of the energy loss depends on the jz-meson energy and the penetrated 
mass and is found, for example, in Heisenberg, Kosmische Strahlung, p. 577. The 


losses are approximately independent of energy in the actual energy range from 0.1 
GeV and upwards with: 


(45) 


The penetrated mass Am between the height of meson production z and a telescope 
at sea level is: 


Am= cai E —exp ( - i)| g/cm? (46) 
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and the loss of 4-meson energy consequently: 


d 
AW,= Wits F exp (—7)| GeV. (47) 


dm cos B 


A telescope records mesons coming from a given height z if the initial kinetic 
meson energy is larger than the total energy loss, which includes losses both in the 
atmosphere and in the recording apparatus. Thus, all mesons of the given energy 
are either recorded by the telescope or not depending on the angle of incidence f. 
This is allowed for by introducing a function J(z,B, Wu xin), defined by: 


J =1 when Wuxi, >AW,+AW, GeV (48a) 
J =0 when Wyyin<AW,+AW, GeV (48 b) 


Wu xin 18 the initial energy of a w-meson starting towards the telescope, A W,, is 
given by eq. 47 and AW, is the energy loss in the recording apparatus (0.1 GeV 
assumed). 


Summary.—The complete expression for the counting rate of a telescope, situated 
at sea level and pointing in a direction making the angle f with the vertical, is given 
by eq. 49a, when the telescope records -mesons which are produced at the height z 
in the interval dz. The telescope has an infinitely small solid angle dw and its infinitely 
small surface da is perpendicular to the direction of dw: 


a ee ghee ; 
dn, =dadal (6, Wo) Ti cos Bp’ exp | oaeTboce fi exp ( alte 
ee] 


R 
Z 
Sea ea ie dw’. 1 
< exp | CW ae dzn dw (49 a) 


D 


I(6,W,) describes the angular intensity distribution of mesons at the production 
level z. R gives the production rate at this level for a primary beam of particles sub- 
tending the solid angle dw’ and coming from a direction defined by the angles «’ 
and f’ (fig. 10). D represents the decrease in meson intensity due to decay, and J 
represents the absorption. The product of the functions [(0,W,), R, J, and D is 
denoted K(6,W,,z) in the following. Thus: 


d‘n, =dadw K(0,W,,z) dz ny dw’. (49b) 


The relation between the measuring direction of the telescope («,/), the direction 
of the primary beam («’f’) and @ is found from fig. 10 and given by: 


cos 9 = sinf sinf’ cos(« — a’) + cos cosf". (50) 


The effect of decay and absorption is that the intensity distribution measured 
deep down in the atmosphere has a sharper maximum than that measured just 
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Fig. 19. The contributions of different production levels to the counting rate of a zenith point- 
ing telescope with infinitely small solid angle dw. The curves correspond to the z7-meson energy 
vn = 3.35 in the C-system, which corresponds to m = 1.00. 


below the level of production. This is due to the energy distribution of the “-mesons. 
Mesons starting in a direction making a large angle 9 with the primary beam have 
smaller energies and will consequently decay more rapidly, seen from the laboratory 
system, than those starting in directions which differ little from the primary direc- 
tion. Furthermore, they are absorbed after a shorter path length than mesons of 
higher energies. 

In fig. 19 we have plotted the counting rate of a zenith pointing telescope with 
the surface da and solid angle dw for different production levels as a function of the 
angle B’ of the incident primary beam. f’ equals 0 in this case. The different curves 
represent the counting rate of the telescope for different production levels. Here it is 
assumed that the production rate is the same for all levels. The uppermost curve is 
proportional to the function J(6,W,) and corresponds to production at sea level. 
With increasing height to the meson producing layer the counting rate is reduced by 
the decay and at a certain angle there is a cut-off when the absorption limit is reached. 
The curves in fig. 19 are calculated for the case m =1.00 which corresponds to a 
a-meson energy of 3 GeV (giving a mean kinetic energy of 2.3 GeV to the u-meson) 
in the forward direction in the L-system (see fig. 18). 

The total counting rate for any f’-value is obtained by adding the contributions 
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Fig. 20. Dashed curves: The counting rate of a zenith pointing telescope as a function of f’, the 
zenith angle of incidence of the primary beam. The curves correspond to three different 7-meson 
energies in the C-system. 

A, The sum of the three curves (i.e. the image when geomagnetic deflection is not taken into 
account). B, The normalized intensity distribution J(#,W,) just below the point of a nucleon- 
nucleon collision. 


from all production layers, i.e. by adding the curves of fig. 19, each of them with a 
weight given by the production rate at the corresponding height and by the angle 
f’, as in fig. 11. 

The same calculations have been repeated for the z-meson energies y,, = 10.00 
and 2.00 and the three distributons have then been added with the same weight, 


according to section 2.1.2. 
The result is shown in fig. 20. We notice that the width of the resulting distribution 


is decreased by the influence of the atmosphere and the decay. 


2.1.3.3. Scattering of u-mesons 


The mesons are scattered by elastic collisions in the atmosphere. The multiple 
scattering is described by the mean square angle of scattering <6”) which is obtained 
by using the Rutherford scattering formulae; cf. Rossi [18, p. 66): 
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where p, =momentum of a particle at the apparatus, p, = momentum of a particle 
at the production level. W; =21 x 10° eV and X, is a constant of the penetrated 
matter which equals 38 g/cm? for air. dp/da is the momentum loss per unit of 
penetrated mass and f is the particle velocity in units of the velocity of light, c. 

We express eq. 51 in terms of energy because the previous calculations have been 
performed in terms of energy. By using equations 52 and 53: 


Cpa Win ae ita VY (52) 
Wein t Wo = me’, (58) 


where W,,;, = kinetic energy of a particle and W, its rest energy, we get: 
We kin 
(62) = AS oe | (W xin + Wo)? 
X (—dpjdx) J [(Wxin+ Wo)? — Wol 


1kin 


d Wein: (54) 


dp/dx has been taken as a constant in the momentum range considered. 
The solution is: 


Wy kin 
W: 1 { Win + W, 1 2Wo+ W. 
02 , kin 0 0 kin | | 
ae STE ME ae oe am | Won )j 8) 


2 kin 


V <0?) has been calculated for different values of Wixin and W,,;, and the result 
plotted in fig. 21 with (—d Win/d x)= 2 MeV/cm. 

Rossi [18, p. 69] has discussed the scattering of a parallel and infinitely thin 
beam of particles, all of the same energy, which enters a scattering substance. 
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Fig. 21. The mean square angle of scattering in air for ~-mesons of different energies. W, ,,,, and 
W, xin are the energies at the beginning and the end of the path respectively. oa 
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Let the direction of the incident beam be along the z-axis in a system of cartesian 
coordinates. The angular distribution of particles, irrespective of lateral displacement, 
is then represented by: 


Q(z, by) = : mat i) 
Ban, | <6?) 
Eq 56 gives the number of particles with directions of motion between #, and #, + dd, 
taken over all values of the lateral displacement y. 3, is the projection of the scattering 
angle in the z — y plane. 

For a parallel beam of infinite extension in both x- and y-directions (which is the 
case in the present calculations) eq. 56 gives the angular distribution at any point 
and is used in the following when estimating the angular spread of the j-meson 
beam. 

This could be done thoroughly by calculating the scattering of mesons which origi- 
nate from different production levels and directions and thus finding the resulting 
angular distribution of “-mesons at sea level. 

However, it is enough for our purpose to show that the influence of the scattering 
can be neglected in the present calculations. 

The angular distribution at each production level, the function J(6,W,), can be 
represented by a Gauss distribution: 


(56) 


I(6, W,) + A exp(— 62/66), 


where 6, is found to be about 16° (Curve B in fig. 20). If we consider the energy case 
y2 = 3.35, the mean kinetic energy of the produced -mesons is given by curve B 
in fig. 18. The contributions from different production levels to the counting rate at 
sea level come from heights up to 20000 m, as seen from fig. 19. The production 
rate increases approximately linearly with the height, see fig. 11. We thus take the 
production level 12000 m as representative of the whole range of production levels. 
The energy loss A W,, in the atmosphere is about 1.7 GeV for ~-mesons coming from 
this height in the vertical direction. From figs. 18 and 21 it is found that V6?) 
varies between 4° and 10° where the high value is at the absorption limit. If we 
take the scattering to be constant over the whole region and equal to 10°, the original 
distribution is widened and the new distribution described by 6; = V 102 + 162 = 19°. 
Thus, the broadening of the distribution is certainly less than 3° in the energy case 
y,, = 3.35. The resulting distribution when adding the three distributions correspond- 
ing to y, = 2.0, 3.35 and 10.00 is influenced still less, as the scattering in the energy 
case y = 10.00 is much less due to the higher “-meson energies. The contribution 
in the energy case y,, = 2.0 is negligible. 

We thus conclude that the influence of the Coulomb scattering can be neglected 
in the present calculations compared with other approximations. 


2.1.3.4. Magnetic deflection of the u-mesons in the atmosphere 


The mesons are deflected in the earth’s local magnetic field between the production 
level and the telescope. The magnitude of the deflection depends on the meson energy, 
the measuring direction of the telescope and the magnetic field at the place of 
observation. At the equator where the magnetic field is about 0.3 gauss the total 
deflection is found to be about 5° for jz-mesons, produced at a height of 20000 m, 
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Fig. 22. Symbols used when caleu- 
lating the magnetic deflection of 
secondary cosmic ray particles. 


z Geographic 
South pole 


arriving at the earth from the zenith direction and having an energy of 1 GeV when 
they reach sea level. As about half of the “-mesons are positively and half of them 
negatively charged [35] this deflection means a spread of the meson flow perpendicu- 
lar to the direction of the magnetic field, and changes the angular intensity distribu- 
tion discussed in 2.1.3.2. 

The deflection of low energy particles in the earth’s magnetic field has been studied 
by Johnson [33]. Calculations of the deflection of secondaries at seal level have been 
made by Bowen [34]. 

In the following, an approximative formulae for the magnetic deflection is 
deduced by taking the particle energy as constant during the motion through the 
atmosphere and equal to the mean value of the energies at the production level and 
the recording apparatus. 

The orbit of a positively charged particle moving towards the earth is calculated 
by treating the problem of a negative particle of the same mass moving in the opposite — 
direction. Fig. 22 shows the situation. The z-axis in the x, y, z system is directed in 
the zenith direction at a given place of observation. The magnetic field B is parallel 
to the z’-axis of a a’, y’, z’ system. The a’- and x-axis coincide and the 2’, y’, 2’ system 
is turned an angle (90° —6) around the x-axis. 6 is the dip angle of the magnetic 
field. The direction of the orbit at the telescope is represented by 7’ and the angles 
a and f define this direction in the x, y, z system. # and y are the corresponding angles 
in the a’, y’, z’ system, and s is the path length taken from origin. As the particle 
energy is assumed constant during the motion, the orbit is a spiral on a cylindrical 


surface. The projection of the orbit in the 2’, y’ plane is a circle with the radius 
of curvature r,: 


r, =7) sin y, (57) 
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Fig. 23. 


where rT) = p/eB =( V/Bc)/V1 + 2mgc?/eV. 


The following relations are found: 


cosy = sinf sina cosd + sind cosf, (58) 
eee sin 6 cos (os B. (59) 
cos 6 sin y 
The particle velocity is v. According to figs. 22 and 23 we have: 

s=vt, (60) 
‘s cos(ka —y) =", 9, (61) 

t 
z=|v cos B; dt, (62) 

0 


where 2 is the particle velocity at the time ¢t and /; is the angle between v, and Z. 
We find from fig. 23 that: 


cosf; = sind cosy — cosd siny sin(? +9), (63) 


which together with eq. 62 yields, as v=»: 
t 
z= | v[sin 6 cos y — cos 6 sin y sin (+ p)] dt, (64) 
0 


where ¢ is given by equations 60 and 61: 


gy =(vt/r,) siny. (65) 
The solution of eq. 64 is: 
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z=1, siny cosd[cos(# + 8/79) — cos] + 8 sind cosy, (66) 


which gives the relation between the height of the production level z and the actual 
path length s. 


For positively charged particles coming from the zenith direction (8 =0) at a 
place where 6 =0 eq. 66 is reduced with equations 58 and 59: 
z=1ry sin(s/ro) ~s8[1 —§ (8/79)")- (67) 
The direction of the orbit at a certain point is given by: 


a =at+Aa, Bb, =B + AB, 


where Ax and Af are due to the magnetic deflection. We find from fig. 22 that: 


cosf, = dz/ds, (68) 
cos a = (da/ds) (1/sinf,). (69) 

Derivation of eq. 66 gives: 
dz/ds= cos, = sind cosy — siny cosd sin(# + 8/79). (70) 


The following expressions for «, is found in an analogous way: 
sin y cos (# + 8/79) 


sin f, 


COS & = 


(71) 


For positively charged particles coming from the zenith direction (f =0) and with 
6 =0 eq. 70 is reduced to: 
cos; = cos (s/79) (72) 


or By = 8/19; 
which is obvious. 

For a “-meson starting at 20 km height and reaching the apparatus with an energy 
of W,=1 GeV the angular deflection is found to be 4.6° using eq. 70, a mean 
energy of 1.9 GeV and a magnetic field of 0.3 gauss. 

Bonnevier has taken the energy dependence of the covered path length into account 
and the more exact formulae are deduced by him and given in appendix. A 
comparison between values of the angular deflection of “-mesons calculated using 
both the approximative and the more exact formula is shown in the following table; 
where the mesons have been assumed to have an energy of 1 GeV when they arrive 
at the telescope: 


Hight of “-meson Angular deflection 
production Approx. Exact 
(m) formula formula 
30 000 7.5° 6.6° 
20 000 5.2° 4.6° 
15 000 3.9° 3.7° 
10 000 2.9° Soke 


The approximate formula is good enough for our purpose when discussing the 
magnetic spread of the “-meson flow at the equator. The deflection as a function 
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Fig. 24. The angular deflection of 4-mesons in a magnetic field of 0.3 gauss perpendicular to the 
motion. The mesons start at different heights and arrive at the telescope from the zenith direc- 
tion. Their energy is constant during the motion, and equal to the mean value of the energies 
at the production level and the telescope. 


of the -meson energy for mesons starting at different heights is plotted in fig. 24. 
The meson energy is assumed constant during the motion. 

The effect of the magnetic deflection upon the image recorded by a zenith point- 
ing, infinitely narrow angle telescope at the equator is estimated when the source 
is assumed to lie in the earth’s equatorial plane. The angular intensity distributions 
of u-mesons orginating from different production levels and recorded by the tele- 
scope as shown in fig. 19, are used together with the corresponding j-meson energies 
(fig. 18). The angular deflection of mesons of different energies between the produc- 
tion level and the telescope is calculated for different production levels and is found 
to vary, for instance: 1°-4° at 10000 m, and the meson energies corresponding to 
the energy case y,, = 10.00. The lower value corresponds to “-mesons which move 
in directions which deviate little from the primary direction (0 0). 

The angular deflection A@ is thus a function of 0, the angle between the primary 
particle direction and the direction of the produced s-mesons: 


Ad =0 (6). 
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particles/sec 


Fig. 25. (a) {(6) is the intensity distribution recorded by the telescope when the earth rotates, 
assuming no magnetic field. A = (0) is the magnetic deflection of the mesons. F'(@) is the 
recorded intensity when the magnetic deflection has been taken into account. 

(b) This shows the path of a u-meson without (left) and with (right) the influence of the magnetic 
field. 


The angular intensity distribution at sea level, corresponding to a certain level 
of production and modified by the deflection of the mesons, is found by the following 
consideration, see fig. 25. 

{(0) is the given distribution and F'(0) is the distribution modified by the magnetic 
deflection. We consider in fig. 25a the element F'(0) 60 which is called C. This element 
represents the number of particles recorded by the telescope during the time interval 
when the earth rotates through the angle 66. Particles within the element C originate 
by deflection from the two elements A and B of the function f(6#). The elements © 
A and B are at the angular distances A 4 and A@z respectively from the element C. 

First we must find the relationship between the shape and position of the element C 
and the shapes and positions of the two elements A and B. Using fig. 26 we deduce 
the relationship between an original element /(#) 66 and the two parts into which 
this element is divided by the magnetic deflection. Fig. 26a shows the original 
element; fig. 26 shows the variation of magnetic deflection with 0; and in fig. 26c 
are shown the two parts into which the original element has been split. It is seen 
from fig. 26 that the angular extensions along 6 of the two resulting parts are, respec- 
tively, 

60{1 —d(A6)/d0} and 66{1 + d(A6)/d6}. 


From the above it is clear that the angular intervals of the elements A and B in 
fig. 25 are 60/(1 —ka) and 66/(1+ kz) respectively, where k = d(A6)/d0. 

Subscripts refer to the derivatives at A and B. As half the number of the par- 
ticles are deflected in both directions from both A and B, we get: 


#(0)0= 4 | 10-009) So + 100+ A6,) | (73) 
A 
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Fig. 26. (a) {(@) is the measured intensity 
at sea level if the magnetic deflection of 
the m-mesons is neglected. (b) The varia- 
tion of the angular deflection with @ for 
mesons coming from a given production 
level. (c) The two elements shown are 
those into which the original element, 
{(0) 60, of figure 26a is divided by the 
magnetic field. 


@ +60 -(69- 88) + no) 


d(48) 
dou +49) 


(In the actual calculations we always have k <1.) 


f= AM) , 16+ 054), (74) 


oe Fo=3| l+ks 1k, 


If {(0) =0 for 0 >6,, where 6,, is a cut-off value for the given function /(#) we get: 


_ HO- As) A 
F(6)=4 i for 0,<0<60n+A0, (75) 
and F(6)=0 for 0>6n+A64. 


Furthermore we have from fig. 25a the conditions: 
0(0 + AO4) =ABa, (76) 
0 (0 —AOz) = ADs. (77) 
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Fig. 27. The change in the angular intensity distribution due to geomagnetic deflection of the 
yu-meson flow. Full drawn curves show undeflected (compare fig. 20) and dashed curves deflected 
distribution. 


The influence of the deflection has then been estimated by introducing the following 
two approximations: 


b(0) = +k6 +1, (78) 
f(0) =a exp(— 6?/b?). (79) ¥ 
This gives: 
wof-B(-8%N) L288) 
a b2 LE b2 1-k 
ro-s| Lok ia Lae (80) 


where A@ is the deflection at 6. 


AJI(6) = f(8) — F (6) has been calculated for each production level and the contri- 


butions from all levels added with the proper weight for the energy cases y, = 3.35 
and 10.0 as described in section 2.1.3.2. 
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The found difference has been added to the intensity distributions obtained in section 
2.1.3.2 and the resulting distributions are shown in fig. 27, where also the original 
distributions have been plotted. 


The influence of the deflection on the resulting distribution, X,’, is not very 
important and has been neglected in the calculations. 


2.2. The polar diagram of telescopes of finite dimensions when the geomagnetic 
influence on the primary particles is neglected 
It is now possible by using the integral eq. 49 of section 2.1.3.2 to calculate the 


polar diagram of a telescope of finite aperture. With the abbreviations in section 
2.1.3.2 the integral is written: 


d‘n, =dadw n) dw’ K (6, W,,z)dz, (81) 


fk(a Ba’ B’ W, z) dz 


°o 


1 


x 


-2 
o ‘da dw-n,dw’ particles/sec. 


“40° §-30° -20° 


Fig. 28a. Polar diagrams of a telescope with infinitely small solid angle, when the geomagnetic 

influence is neglected i.e. the intensity measured by a telescope da dw tilted at different angles f to 

the vertical as a function of f’, the zenith angle of the primary beam. The azimuth angle of the 
telescope, « = 90°. 
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Fig. 286. As in fig 28a, with a = 67.5°. 


where @ is a function of «, B, «’ and f’ according to eq. 50 and dw =sinB da df. 
For a zenith pointing, cylindrical telescope at sea level the counting rate due 
to the primary beam is: 


22 “max oo 


Any, = Ny) do’ [ i} sin Bb dadB [da cos B { K (6, W », 2) dz. (82) 
=0 B=0 a 270 


a= 


The effective measuring surface is (compare 1.4.1.): 
[cosp da =a cosf G(f). 


The rest of the integral has been solved with numerical calculations. The integral 
was first solved for different measuring directions in zenith as well as in azimuth for 
a telescope with infinitely small solid angle and for different directions of the primary 
beam. The result thus obtained is the polar diagram of a telescope with infinitely 
small solid angle pointing in different directions and situated at the equator, the 
source being in the equatorial plane. This is shown in figures 28a and 286 for two 
azimuth directions and for a number of zenith directions. The polar diagram of a 
finite telescope has been obtained with the help of these and similar diagrams. 
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Fig. 29. The contributions to the counting rate of different production levels z and two different 
7-meson energies in the C-system. 


The integration giving the polar diagrams in fig. 28 has been done as follows: The 
atmosphere was divided into a number of horizontal layers of finite thickness Az = 
2000 m. For each layer a set of curves was calculated which give the counting rate 
of the telescope pointing in different directions when recording w-mesons which have 
started towards the telescope with different energies from the center of the layer. 
For a certain z-meson energy in the C-system these energies are given by the angle 0 
between the primary particle beam and the direction of the meson flow produced as 
discussed in section 2.1.2. As an example a few diagrams corresponding to 7-meson 
energies y,, = 10.0 and 3.35 in C-system are shown in fig. 29. Both the decay and the 
absorption have been taken into account (the functions D and J of eq. 49). The 
production rate represented by the function R is however still put equal to one. 
Diagrams of this type have been computed for three 7-meson energies in the C-system, 
namely y,, = 10.0, 3.35 and 2.00. As these energies were chosen so as to represent the 
same number of mesons within the corresponding energy intervals (compare 2.1.2) 
the curves in fig. 29 belonging to the same atmospheric layer have been added with 
the same weight. The combined diagram represents, approximately, the total count- 
ing rate of the infinitely narrow angle telescope due to meson production within 


the given layer. 
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Fig. 30. The counting rate 
of zenith pointing telescopes 
with different apertures as 
a function of the zenith angle 
of incidence f’ of the primary 
beam. As f’ equals Y’, the 
longitudinal position of the 
source, these are the polar 
diagrams for D = 0 when the 
geomagnetic influence is 
neglected. These functions 
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The next step was to summate the contributions from all layers in the atmosphere 
with regard to the variation of production rate. For a given direction of the primary 
beam («’, 6’) and measuring direction of the telescope («, /) the angle 6 was obtained 
from eq. 50. The contribution from each layer was then calculated by taking the 
counting rate from the combined diagrams of fig. 29, multiplied by the value of the 
function F for the corresponding layer. This was repeated for all layers up to a height 
above which there is no contribution to the counting rate. 

The procedure was repeated for a number of measuring directions « and f and 
different angles of incidence f’ of the primary beam. As the telescope was assumed, 
for simplicity, to be situated at the equator, «’ of the primary beam is always 90° 
and f’ equals the longitude angle ‘’ of the source. We thus finally arrive at the 
diagrams of the type already shown in figures 28a and 286, which have been com- 
puted for B =0°, 10° ... 80° and «= +90°, +67.5°, +45.0°, +22.5°, 0°. These 
curves are solutions to the integral of eq. 82: 


[K (0, W,,2)dz. 
2=0 
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Fig. 31. The percentage increase of the counting rate when measuring a very small source element 
with telescopes of different apertures. The cos? f law is assumed to hold at sea level for the intensity 
distribution of the “background”’. 


The case of a finite telescope can now be tackled using these diagrams and requires 
an integration over da dw in eq. 81. Our cylindrical telescope was assumed to point 
in the zenith direction. Its total solid angle was divided into a number of sections 
both in zenith and azimuth, each section being centred at those values of « and f 
for which the counting rate of an infinitely narrow angle telescope had been computed. 
The total counting rate for a given direction f’ of the primary beam was then found 
by summing up the contributions from all the parts of the solid angle. The values of 
the counting rate for each part were taken from the above mentioned diagrams and 
multiplied by the value of the function cos u G(u) for the mean direction f of the 
section in order to allow for the effective measuring telescope surface in the same 
way as described in section 1.4. The result is plotted in fig. 30. 

When applying these results to the measurements of daily variations it must be 
pointed out that what we actually observe, when measuring with cosmic ray tele- 
scopes, is the ratio between source and background intensities and not the source 
intensity itself [36]. If we were to measure an amplitude, of say 1%, in the daily 
variation with an infinitely narrow angle telescope, assuming the source which 
causes this daily variation to subtend an infinitely small solid angle, then a broad 
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angle telescope measures a smaller amplitude, because it records a higher ‘“‘back- 
ground”. hed 

According to fig. 28a a telescope (da dw) pointing in the zenith direction records a 
maximum counting rate of 


(9.45/m)dadaw ny Aw’ particles/sec 


due to the source which subtends the small solid angle Aw’. If we assume the ratio 
between the counting rates of source and background intensities to be 1% we get: 


ny Aw'/I, = 2/945. (83) 


The counting rate of cylindrical telescopes due to the background is calculated in 
appendix IV assuming the cos?f law to hold. The result is plotted in fig. AIV.1. 

With the relation in eq. 83 and with fig. AIV.1. the relative counting rate of source 
and back ground intensities recorded by telescopes of different apertures has been 
calculated from fig. 30. The result is shown in fig. 31, where the big difference between 
a narrow angle telescope and an ionization chamber should be noticed.(') 


2.3. The zenith angle dependence of the radiation at sea level 


The variation of the counting rate of a telescope with the zenith angle f is found 
experimentally to be approximately a cos?/ curve at sea level [38 and 39]. This 
fact offers a valuable check of the calculations as regards the atmospheric part 
and has been carried out for a narrow angle telescope assuming isotropic distribution 
of point sources over the hemisphere using the described calculations. The result is 
shown in fig. 32 where the cosf, cos?f and cos?f curves have been drawn for compari- 
son. The calculated curve agrees fairly well with the cos?f curve. The rather good 
agreement indicates that the assumptions regarding the meson production and the 
primary energy of 20 GeV give an essentially representative picture of the phenomena 
in the atmosphere. 
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(1) Amplitudes of the order of 1 % in the daily intensity variation using narrow angle telescope 
have been reported by Sarabhai [37]. Such results could be explained by these calculations. 
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A few remarks should be made in connection with the zenith angle dependence 
of the radiation. The first one applies to ionization chambers. It is well known that 
the maximum contribution to the counting rate of an ionization chamber per unit 
of the zenith arc comes from directions between two concentric cones centered around 
the zenith direction and defined by the apex angles 28 and 2(8 + dB) where B = 35° 
if the cos?# law holds for the radiation. It is obvious that this direction has no prac- 
tical significance when discussing the directional properties of the chamber. Its most 
sensitive measuring direction is in the zenith direction and only in a hypothetical 
case when the whole particle flow is varying with time in the same way in all direc- 
tions can we talk about f = 35° as the most sensitive direction. 

Another remarkable point is the following. A narrow angle telescope is tilted at 
an angle f to the zenith direction. Its counting rate due to the background radiation is: 


y = N, cos?f  particles/sec, 


where J, is the counting rate when pointing in the zenith direction. A change in 
means a change in the measured intensity according to: 


Av/y = — 2 tg B AB. 


A percentage change in the counting rate of say 0.1% at 6 =45° corresponds to a 
change in Af of 0.03°. Daily changes in the measuring direction of this order of 
magnitude could introduce considerable errors in the image of a cosmic ray source. 
It is consequently important that a tilted telescope is mounted in a stable way. 


3. THE COMPLETE POLAR DIAGRAM OF A COSMIC RAY TELESCOPE 


At this stage we are able to discuss the complete polar diagram of a cosmic ray 
telescope taking into account both the geomagnetic deflection and the influence of 
the atmosphere. 

The angular position of a source element on the sphere of asymptotic directions 
is given by ’ and Q; the position determines the direction of incidence («’,6’) of 
the primary beam at the top of atmosphere (see fig. 10). The intensity of the 


beam is 
dn,=I(p)dp_ particles/em*-sec:sterad, 


and the beam subtends the solid angle (cf. eq. 1I8b) Aw’ = M-AQ where AQ is the 
solid angle subtended by the source element. /(p) represents, as before, the primary 
differential momentum spectrum of the source and M the magnification of the 
geomagnetic field. 

Secondaries are recorded by the given telescope provided the source element is 
within the projected extended solid angle of the telescope on the sphere of asymp- 
totie directions for a certain momentum p. This solid angle extends beyond the 
projected total solid angle defined in section 1.4 due to the influence of atmosphere. 
The extension is found from fig. 30 when considering primary particles of momen- 
tum 20 GeV/c. 

The beam intensity was introduced in section 2.1.1, and the counting rate of the 
telescope at sea level is found from eq. 81 in section 2.2 when the atmospheric in- 
fluence is taken into account. The contribution to the counting rate from a surface 
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element da of the telescope and solid angle dw due to “-mesons produced by primaries 
within the momentum interval dp is consequently: 


dn, =I1(p)dp Aw’ da cos ud | K(6, W ,, 2) dz. (84) 
0 


The direction of dw is given by « and f. As the cylindrical telescope points in the 
zenith direction we have u =f. 6 is a function of «, B, «’ and f’ given by eq. 50. 
Integration over the whole telescope surface yields: 


d?n,=I(p) dp Aw’ acosB G(B) dw { K (6, W ,, 2) dz, (85) 
6 


which is the counting rate of the telescope for ~-mesons arriving within the solid 
angle dw (compare eq. 20). 
Integration over the total solid angle of the telescope gives finally: 


22 “max 0 
dn,=Tp)dp MAQa| f | cos B O(p) si Bdaxdp f K(0, Wy.) a2} (86) 
0 


a=0 B=0 


The function within the brackets, Q(B’, Wp, Umax), Was solved for telescopes with 
different apertures Umax in section 2.2 and plotted in fig. 30. The counting rate 
due to primaries in a momentum interval dp at p is obtained from eq. 86. 
The value of f’ is given by the angular position Y’ and ® of the source element on 
the sphere of asymptotic directions. (The value of «’ is of no importance as the 
telescope is assumed to point in the zenith direction.) M is taken from diagrams of 
the type shown in fig. 8b. Contributions from all momentum intervals for which 
Q(f’, Wy, Umax) #9 give the total counting rate and thus the value of P at ¥’, ®, 
giving the complete polar diagram wken both the geomagnetic deflection of the 
primaries and the influence of.the atmosphere have been taken into account. 

Let us however discuss eq. 86 once more with particular reference to the function 
K(0,W,,z). As MAQ is independent of both « and f, eq. 86 is also written: 


22 “max 


dn,= | | acosBG(f)sinBdadB I(p)dpM AQ { K(a,B, a’, fp’, W,,2)dz, (87) 
0 


2-0 p=0 


where I(p)dp MAQ{ K(a,B, a’, B’, W,,z) dz (88) 
0 


represents the number of «-mesons arriving at the telescope from the direction «, p 
produced by primary particles which come from the direction «’,f’. The integral 
was solved in section 2.2 for W,=20GeV and plotted in fig. 28. It relates the 
number of “-mesons measured in direction «,6 at sea level to the number of 
primaries at the top of the atmosphere coming from the direction «’,p’. 

A function which relates the measured number of one sort of secondary particle, 
for example jz-mesons to the number of primaries within different momentum inter- 
vals, is called “the over-all multiplicity’ by Neher [40] or ‘yield function” by Simp- 
son et al. [41]. 

Assuming the earth’s field to be a dipole field, and by using this field as a magnetic 
analyzer for the arriving primary particles, it has been possible to deduce the yield 
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Fig. 33. (a) The contributions of different primary momentum intervals to the polar diagram 
along ® = 0 of an ionization chamber at sea level in Huancayo. (b) The polar diagrams along ® = 0 
for teleseopes with different apertures at sea level in Huancayo. 


function from measurements of the latitude effect. This has been done by Neher 
[40] for the ionizing component and by Simpson e¢ al. [41] for the nucleonic compo- 
nent; the yield functions obtained are those of the isotropic radiation. The yield 
function Y(p,f) for the ionizing component derived by these authors is related to 
the one given in eq. 88 by: 


2x n/2 oo 
1 
Y (p, B)= aE: | | sin B’ da’ ap’ | K (a8, OC Bis We) dz: (89) 
ee p’-0 6 
Y (p,B) is a mean function taken over all primary directions «’ and /’ in the hemis- 


phere which are not forbidden for particles with momentum p. 
The yield function derived for an ionization chamber is still more complex: 


2a - n/2 
1 i 
Fi) = 5 | i Y(p, B) sinBdadp (90) 
a=0 £B=0 
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The yield function times the primary spectrum is called the effective primary 
spectrum. The effective primary spectrum of an ionization chamber defined in this 
way is known up to 13 GeV/c. Fonger [42] has deduced a spectrum above this point 
of the form (1 + W,) 1”. The correct way of calculating the polar diagram is, how- 
ever, to use eq. 86 taking into account the variation of W, in the function K (0, W,,2). 

This energy dependence means roughly that the energy available for meson produc- 
tion increases with the primary energy and also that the meson energies in the 
forward direction in the laboratory system are increased. Thus the mean energy of 
the z-meson flow is an increasing function of the primary energy. Due to absorption 
and decay, fewer jz-mesons reach sea level per given number of primaries in the lower 
primary energy range. This corresponds to a decrease of the yield function with 
decreasing energy. 

The variation of W, has been taken into account approximately by using the 
Fonger-spectrum (1 + W,)*” when calculating the complete polar diagram at sea 
level shown in fig. 33. The error introduced by this approximation is probably of 
minor importance at present compared to the errors introduced by our lack of knowl- 
edge about the geomagnetic deflection in the low momentum range. It is also 
assumed that the function Q is independent of W, and the solutions for W, = 20 GeV 
have been used. This assumption introduces an error which is of little importance 
for the shape of the polar diagrams when considering apertures larger than 15° and 
energies above 20 GeV. 

It is, however, possible to solve the equations presented in this paper for other 
values of the primary energy than 20 GeV and thus take the energy dependence 
into consideration in a more correct way than by using the effective primary spectrum 
integrated over all primary directions. 

The polar diagrams of fig. 33 are of fundamental importance when calculating 
the image of a cosmic ray source lying in the equatorial plane. On the other hand 
when the image is given, the polar diagrams are used to correct for the influence of 
the geomagnetic field and the atmospheric spread, thus finding the real angular 


position in space and the shape of the source. This is further discussed in section 3.1 
and chapter 4. 


3.1. The image of a source element obtained by using the polar diagram 


Consider a source element infinitely far away from the earth in the equatorial 
plane. Its angular position is given by ” which is an easterly longitudinal angle 
reckoned from the sun-earth line, see fig. 346. When the earth rotates the image — 
of the source is given by the polar diagram P(‘Y’,® = 0) as in fig. 34a. Since ¥ is 
taken as increasing in an easterly direction, we find that the shape of the image is 
the reverse of that of the polar diagram. 


4. THE RESOLVING POWER OF COSMIC RAY TELESCOPES 


As in ordinary optics, the resolving power of a cosmic ray telescope could be defined 
as the smallest distinguishable angular distance 6(wt) between two cosmic ray source 
elements, as in fig. 35. The criterion of observability could be the statistical signifi- 
cance of 6 N/N, in the figure. This should mean that the resolving power depends on 
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d [nm (‘£")] = km: d QP (4) 


d [nm (@ t)] =km-d Q cos(me) P(Fo- atte) 
particles/sec 


Fig. 34. The image of a source element 
recorded by a telescope using its polar 
diagram. The fulldrawn curve is the 
image of a source in the direction ¥”. 


12 0 
The dashed curve corresponds to a 
source in the direction ¥’ +¢. dQ is 
the solid angle subtended by the VV 
source element. km is proportional to 18 
the amplitude of the mth intensity € 
distribution and to the surface of the 


telese ; . 
Saas Im Im:cos (mé) particles/cm* sec: sterad 


the geomagnetic position and aperture of the telescope which together determine 
the polar diagram and consequently the image of a source element. It also depends 
on the size of the telescope measuring surface as this, together with the aperture, 
determines the magnitude of the counting rate. For a given surface, a wide aperture 
gives a high counting rate, and the statistical determination of 6N/N, is more ac- 
curate than with a narrow aperture and the same surface. But a wide aperture means 
also a polar diagram with large extension along the ‘’-axis as has been shown in 
the preceding chapters. The image obtained is consequently broad which reduces 
6N/N,. Thus a compromise is required; the optimum aperture being found empiri- 
cally. 

A more practicle definition of resolving power is the capability of the telescope to 


Fig. 35. The resulting image, x, of two source 

elements at an angular distance 6(wt) apart. 

@ is the angular velocity of the earth. N, is 

the counting rate due to the isotropie back- 
ground, 
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reproduce higher harmonics in the extra-terrestrial intensity distribution. This defini- 
tion is particularly convenient because the measured images are usually analyzed 
in terms of first and second harmonics which describe the angular extension of the 
source around the earth’s axis. 

Let us assume that the particle intensity around the earth varies sinusoidally with 
the period m 2z along the equator with maximum intensity J, in a direction given 
by ¥” which is an easterly longitudinal angle reckoned from the sun-earth line as in 
fig. 34. The longitudinal variation of intensity is represented by a distribution of 
source elements along the equator, all with the same momentum spectrum but with 
a source strength which varies sinusoidally with maximum intensity in the direction 
’. A source element subtending the solid angle dQ produces an image as described in 
section 3.1. The images produced by different source elements are displaced relative 
to each other and all subtend the same angular interval along the time axis wt 
but their magnitude depend on the corresponding source strength. The resulting 
image is obtained by integrating the contributions from all directions «. 

This is done by developing the function P(‘’,® =0) in a Fourier series. Gener- 
ally the image due to the source element subtending the solid angle dQ at the angle 
(Poe) is: 
km dQ. cos(me) P(V5 — wt + €) 

oo 
= k,, dQ cos (me) dP, cos [vy (Wo—wt+e)—W,], (91) 
where Yo =‘ +2, see fig. 34. km is proportional to J,, and to the surface of the 
telescope. 

m means the mth harmonic of the intensity distribution along the equator and P, 

denotes the different Fourier-coefficients of P(Y,® = 0). 


Keeping wt constant we integrate the contributions from all directions ¢ with 
dQ) = de d® where d9@ is the angular extension of the source element along a meridian: 


d[nn(wt)| =k, dD J cos (me) > P, cos [vy (Wo —wt+e)—W,]de. (92) 
0 — 
For the yvth term we have: 


d[rn,(wt)|=k,d® i P, cos (me) cos [vy (¥5 — wt — WY, /v) + velde. (93) 
0 
Integration yields: 


22 


d [my (wt)]=k,d® P, Jeos, | E (m—v)e  sin(m+yr)e 
| 2 2(m—v) 2(v+m) 
. 7 cos(y+m)e cos(y—m)e 
=a 0D, 27 
ro /| 2(v+m) 2(v—m) || ed 
with 0,<» (Ys ot-2). 
, v 


The integral is different from zero only when y= m 


d[n,(wt)|=k,d® P, x cos lm (vi —wt— 1) . (95) 
v 
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Thus the image of a sinusoidal intensity distribution along the equator is a sine 
wave with the same period as the period of the given distribution and with maximum 
amplitude at 


awl 
wot= Fo— > F. (96) 


Furthermore it is found from eq. 95 that the amplitude of the vth harmonic of 
the image is proportional to the amplitude of the yth harmonic of the polar diagram 
P(’,®) when developed in a Fourier series. 

The resolving power is consequently determined by the Fourier coefficients of the 
polar diagram of the telescope. 

The counting rate due to a constant intensity along the equator which corresponds 
to m= 0 is: 

d[ny(wt)] = ky d® P, 22. (97) 


Let us take an example and discuss what an ionization chamber at Huancayo 
measures. The polar diagram of fig. 336 in chapter 3 has been analyzed and the 
amplitudes of the zero, first, second and third harmonics are respectively in the units 
of fig. 33: Py) = 12.5, P, = 17.9, P, =7.9, and P, = 1.3. The maxima occur at Y’, = 54°, 
VY, =2 (55° +n X 180°), and ¥, =3 (38° +7 xX 120°). 

If the intensity distribution along the equator varies as 


I 
> — cos (me), 
m £0 
where > J,, cos(me) represents the anisotropy and J, the isotropic background, 


m 
the relative amplitudes of the measured first, second and third harmonics, found 
form eqs. 95 and 97, are: 


bat UGE Te ie le a a 
sulaly 9 2 Py’ 0 2 Po 
With the values of P, in the actual case we get: 
ey. 0.72, Ly 0.32, ds 0.05. 
ky ko ky 


Thus the amplitudes of the extra terrestrial intensity distribution are reduced approxi- 
mately by these factors when measured with an ionization chamber at Huancayo. 
The approximation depends on the incomplete polar diagram which was used. As 
pointed out in section 1.4.4 its shape is not determined in the lower momentum range. 
The complete polar diagram will reduce the above figures and also change the angles 
of maxima ’,, Y, and Y's, giving an increase in the case of ‘f’. 

In the discussion, we have only considered an intensity distribution along the equa- 
tor but we may also have a sinusoidal intensity distribution along a meridian. When 
the period of the distribution is 27 both along the equator and along a meridian, 
the total distribution around the earth corresponds to a case when the earth moves 
relative to the isotropic distribution of particles. The shape of the image is in this 
case not strictly the same as with a distribution only along the equator and has to 
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be integrated over all values of ® using the corresponding polar diagrams P(‘f’,®). 
But, as the projected solid angles on the sphere are mainly symmetrical around the 


equator and the main contribution to the counting rate comes from this region, the 
shape of the image will not differ greatly from the one treated above and in any 


case its extension along wt will be the same. 


APPENDIX 
AI. Calculation of the effective measuring surface of a cylindrical telescope 


The telescope configuration is shown in fig. AI.1. The shaded surface is: 


S=9r | Pane — 222 Pee (ALI) 
L ; d 
with @max given by: COS Mmax = or tg u. 
We introduce: 
ee sin 29 
G Ne aha a max] | 
(uw) Ae = 9 | (A 1.2) 


Fig. AI.1. The cylindrical telescope considered. 


248 


ARKIV FOR FYSIK. Bd 14 nr 16 


u 


Fig. AI.2. The relative counting rates of telescopes with different apertures as a function of the 
angle of incidence. 


The product cos wu G(u) is plotted in fig. AI.2 for different values of Umax and is 
the relative counting rate of a telescope which makes the angle w with the parallel 
beam of incident particles. 


AIL. The total production rate of secondaries in the atmosphere 


Eq. 32 in section 2.1.1 is written: 


$a 
dn C 
—= — aa +} ‘dp’ ATI.1 
a4 | exp|- (+25) | tee" ap (AIL) 
p*-0 
—2/l 
2 ig ae No M% pen yo ey = 
with A=2nx IL? B=z2/l, C L 
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The integration is carried out for f’ varying between 0 and 47 at a constant height z- 
Thus: 


n/2 


dn _ _B Huet? dp AIL.2 
ia Ae [ exe( ap) 8b dp’. ( .2) 
0 


By the substitution y = C/(cosf’) we get the integral: 
~ —y 
dt aoe | ae 8 (AIL.3) 
y 
Cc 
which is tabulated in Jahnke-Emde [43]. dn/dz is plotted in fig. 13. 


ATI. Intensity and energy distributions in the laboratory system of ~-mesons 
produced in nucleon-nucleon collisions 


The intensity distribution F'(0,W,) in the laboratory system (L-system) of sec- 
condaries, which are emitted at an angle 0 to the direction of an arriving primary 
particle with energy W, is given by Bradt, Kaplan and Peters [30]. Thus 


dN = F(6,W,) d0 (AIIL.1) 


is the total number of secondaries emitted in directions between two concentric cones 
defined by the apex angles 20 and 2(9+d6). 

The number of secondaries emitted within a beam subtending the infinitely small 
solid angle dw is: 


Be CA a 
dn= Se ay dw=I1(0,W,)dw. (A III.2) 
The following symbols are used: 
; sin 6’ 
X=Yp tg 0 ee ate (A III.3) 


The prime denotes quantities in the centre of mass system (O-system). The subscript 
p refers to the proton. 


y=1/V1—B? with B=oJc, (A III.4) 
a2, (AIIL5) 
Wot = 2myc? y, = total energy in C-system, (A IIT.6) 
Wor = Mpc? (1+ yp) = total energy in L-system, (AIII.7) 


Up _ velocity of proton in C-system 


m=- : : ‘ 
v, velocity of z-meson in C-system 


(A IIL.8) 
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e 
‘ 


Fig. AIII.1. The relation between the velocities of a proton 
and a z-meson in the centre of mass system of two colliding 
protons. 


The average energy of mesons in L-system <y,> is related to the average energy in 
O-system <y,> by: bebe 
Ya> = “Yu? Yo- (AIII.9) 


We have to distinguish between three cases: m > 1, m = 1 or m <1 (see fig. AIII.1). 


m>1 
lO.W Nikitin Hafacceipt Vp m2 + 1— (m2 — 1) o&+2mV1 — (m?— 1) 2 
U, ‘per (a2+ 1)? cos?6 V1—(m?—1) a2 
(A IIT.10) 
hast fageeets oe? eee 1 ped Van 
yg = Val¥o) yo + 0 )t Vy WV yp IVI (PVP ery) 
a ue 
The plus sign applies when 6’ <9),. 
The minus sign applies when $2 >6’ > 6),. 
m=1 
N at Yo 
= : _—— AIII.12 
es azsin@ («2+1)? cos?@’ ( ) 
Che N gays yer Vee] 
_ (aly) (Yp + a) + V5 1 Vy 1. (A 111.13) 


ve a+ 1 
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m<1l, 
0<i2 
Te. W eee «yb [m+ ViR= (mt 1ot] (A III.14) 
»?!~4asind o&+1 cos? V1—(m2—1) 
(yay) (ye + a2) + Ve =I Vy = LVI (mV) = a Try 45) 
A aay eke opal! 
06> 32 ae 
f =. oa a 2 
TU Vee ee [m— Vi = (m?= I) oA] (AIIL.16) 
4xsinO «+1 cos?6 —/1—(m?—-1)c2 
_(Walyo) (vot + 02) — Vy = 1V yp LVR (mY? ay 


* +1 
From eq. AIII.3, if m=1, the relation between 6 and 6” is: 


u 6’ 6’ 1/2m,c? 
tg6= 189 8 OV Wea’ 


where Wi; is given by eq. AIII.7. 
Half the number of particles in the L-system are within a cone defined by the 
apex angle 26,. The value of 0, is found from eq. AIII.18 by putting 0’ = 32: 


tg 0, = V2m,c2/Wtot- (A IIT.19) 


A IV. The counting rate of cylindrical telescopes assuming a cos’ £ distribution 
of the intensity 


We consider fig AI.1 and find the expression for the counting rate of a surface 
element da within the shaded surface. The beam of particles subtends the solid angle 
dw and its direction of incidence makes the angle u to the vertical. As the telescope 
points in the zenith direction wu =f. The beam intensity is given by: 


I =I, cos*f particles/em?- sec: sterad. AIV.1 


Thus the contribution to the counting rate from the element da due to particles coming 
within the solid angle dw from the zenith direction f is: 


d®y = I, cos*B da cosB dw __ particles/sec. AIV.2 


This applies to any surface element of the shaded surface. The counting rate of the 
telescope is consequently: 


d?y = I, cos?B Sdw _particles/sec, AIV.3 


where S is given by eq. AI.1 in appendix I. 
Integration over the total solid angle yields: 
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Fig. AIV.1. The variation of counting rate with aperture of a cylindrical telescope assuming the 
intensity to vary as cos?f where f is the zenith angle. J, is the intensity in the zenith direction. 
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yal [ i 0088 (Pax aE) sin dad, (A IV.4) 


a—0 p=0 
where COSYmax =Ctgu and C=d/2r=1/u. 
As « is independent of f we get: 


Pmax 
y=I,4a | cos? 8 (nas Sates) sin B dp. (AIV.5) 
6 


After change of variable (tg 6 =: cos m) and dropping the subscript we have: 


n/2 
a 2 cos Pp : =n *?) ra AIV6 
v tote | cod sin /~p (r =r Q: ( .6) 
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Integration by parts gives: 


n/2 n/2 


4 | 
y=I,4a | et) (»->3 ") I, 4a [ wef) 2sintg dg, 
! 


0 


COS ~ Sin /P 1 1 : 
1 + p? cos? ~)8 ? 4u? (1+ p? cos? p)? 


where f(g) = {, 


Eq. AIV.7 is now written: 


m/2 
IU 2 sin? w 
die eh [ eer 


0 


The solution is (cf. Meyer zur Capellen, Integraltafeln, p. 156): 
vy =a1, 42(1 — costa). 


Eq. AIV.10 has been plotted in fig. AIV.1. 


(A IV.7) 


(A IV.8) 


(AIV.9) 


AIV.10 


Another telescope configuration has been treated by Newell e¢ al. [44]. They have 
calculated the counting rate of a telescope consisting of two parallel Geiger—Miiller 
tubes. Witmer and Pomerantz [45] have discussed the counting rate of a telescope 


with rectangular surface. 
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Magnetic deflection of secondary cosmic ray particles in the atmosphere 
By ByOrn BonNEVIER 


A secondary cosmic ray particle which is formed in the atmosphere moves towards 
the earth and suffers energy loss during its motion. If it is charged, its motion also 
depends on the earth’s magnetic field. The aim of this paper is to calculate the length 
and the angular deflection of a particle orbit. Instead of calculating the orbit of a 
particle moving towards the earth one can obtain the same result by calculating the 
orbit of a particle with opposite charge moving in the opposite direction. 

The problem is as follows: We assume a point at sea level where there is a magnetic 
field B directed at an angle (90° — 6) to the zenith direction. 6 is the dip angle. A 
particle starts in the direction 7 at an angle f to the zenith direction and at an angle 
y to the magnetic field as shown in figure 1. As we treat the case of a particle moving 
upwards, it gains energy during its motion. We want a relation between the length of 
the particle orbit and the height above sea level. The magnetic field is assumed to 
be constant throughout the length of the orbit. Scattering in the atmosphere is 
neglected. 

A particle orbit in three dimensional space can be described by the vector 
|E(s), n(s), €(s)|. The the radius of curvature r is given by 


1 \/ d2é\2 d?n\2 d20\2 
= —2 ] 
r (Fe) + (7 a ds?} ’ (1) 
where s = the length of the arc. 


If the particle moves in a homogeneous magnetic field and there is also a force 
directed tangentially to the particle orbit we have: 


| vy /Y, = k = ctgd; 


vj, = the particle velocity parallel to the magnetic field 

v, =the particle velocity perpendicular to the magnetic field and positive in 
anticlockwise direction 

k =constant 


If the magnetic field is along the ¢-axis we have: 


GS FE 
d&*+ dn? 
and ds? = (1+ k?) (d&*+ dy]. 
dy 1 dé\2? df ek 9 
Thus we get: af 7! Tae (3) Srehpert (2) 


and after derivation and substituting in [1] we arrive at: 


1 dé 1 


r —ds®" 1—(1+ k) (dé/ds)> 
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Fig. 1. Coordinates used in this deduction. 


In order to get 7, the radius of curvature of the projection of the orbit in the 7-plane 
we put: 


1/r, =V d?é/do? + d?n/do?, 
do? = d& + dy? = ds?/(1 +k). 


The derivatives d?&/do* and d?y/do* can be expressed in terms of d?é/ds? and 
d&/ds. We then get: 


1/r, = (1 + ®)/r.(2) (4) 


In order to find another expression for r, we make the following approach. Let 
W(s) denote the kinetic energy of a particle and W, its rest energy. Then the 
particle momentum is: p =[W(W +2W,)]*/c. 


As p/p, =k we get p, =p/V1+k® and 


B.S AWW Ee Wa 


"1 TZelB cl/Ze|B Vig de (8) 


We have now obtained two expressions for r,, (4) and (5). Putting these equal 
we get: 


1 c|Ze|B —. dzé 1+? 
= 1+h?= + - 
rt, ([W(W+2W,)}* V+ ~ 08 V ee (1+ k) (dé/ds)2 (6) 


(*) This result could also have been obtained with the relation between the radius of curvature 
for a surface in a certain direction and the main radii of curvature. 
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The ¢-axis has already been chosen so that it has the same direction as the magnetic 
field. Now we choose the &-axis in such a way that the ££ -plane is a tangential plane 
to the particle orbit at sea level. The origin in our coordinate system is this point 
of contact. For a particle starting from the origin we take ds>0. At the origin 
dé/ds >0 gives the direction of the é-axis and the choice of a righthanded system 
gives the 7-axis. The special case d&/ds =0 at the origin is trivial as the motion 
is then parallel to the magnetic field. We thus choose the minus sign in eq. (6) 
In (6) we make the substitution: 


t= (dé/ds)V1+k 
and get the integrals 


t 


Frag erect =| Ze) B fe (7) 
Vl-# /W(W+2W,) 


0 


where the energy W is a function of the path length s. 

In order to find a relation between W and s, when we treat the case of a secondary 
cosmic ray particle moving through the atmosphere, we consider the motion along 
the short path length ds. The energy gain is dW/ds=o(dW/ods) =oq, where 
o=the density of the penetrated matter. The atmosphere is assumed to 
be in static and thermal equilibrium thus 9 =e,e "", where h = height above sea 
level, /=scale height and 9, =air density at sea level. dW/ods =q is given for 
different energies and particles for example in W. Heisenberg, Kosmische Strahlung, 
Berlin 1953. To simplify the calculation, g is taken as a constant which is a good 
approximation in the relevant energy range. Thus we have 


dW=qoe "ds 
| (8) 


or W=W,+ o,{e "ds 
0 


which enters into the integral (7). 

To solve this integra! we have to make a zero order approximation for the desired 
function h(s). We put h =scosf. It is then possible to integrate and obtain d&/ds 
as a function of s. Using dé/ds we get dh/ds and one more integration gives a first 
order approximation to h(s). In principle we could go further to higher order approxi- 
mations, but the first order approximation is good enough. Now we have 


s Ws 
ds cl/Ze/B ALN 
te Ze) B | ae 
arc cost = ¢/Ze/ aes cosp J VW(W+2Wo)(W3— W) 
0 1 


with the substitution: 


s cos B 
W=W,+(W,- W;) [1 —exp ( ~ ] )|. 
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where W, = particle energy at sea level 
W. = particle energy at path length s 


W, = W, + ¢0ol/(cos) 
W, 22 W, athe pee loss of a particle moving from the top of the atmosphere 


to sea level. 
Integration gives: 


de anaes eS) GUE ee [x cosh = 
ds V1+k  \cosBVW,(W3+2W,) 0 


W,(W,+2W,) ( WatW, Tong }|| 
—ar cosh - 
W (Ws ral W,) die td Dae Wo Wais W,) 


= d cos 9 = sin y cos 8, (9) 
Vl+k 


where @ is the expression within the brackets. 
Now we want to express dh/ds in terms of d&/ds. In addition to the angles Byd 
earlier introduced we need the angles: 


€ =the angle between the ¢ —h and ¢ — & planes; 
gy =the angle between the € —h and h — T planes. 


From equations (2) and (9) 
dé/ds=sin y cos 8, 
dn/ds= +sin y sin 0, 
dl/ds=cosy. 
The plus sign applies for a positive particle moving towards the earth. Since 
h = [cosd cose; cosdsine; sind] 
we get 
dh/ds = (d&/ds) cosdcose + (dyn/ds)cosdsine + (dé/ds)sin 0. 
Using the spherical sine and cosine theorems: 
cosB = cosdsiny cose + sind cosy, 
cosy = sinf cosdcosp + cosfsind, 
sinysiné = sinfsing; 
to eliminate ¢ and y we get: 
dh/ds = (cosB cos?6 — sinB sind cosdsin a) cos6 + 
+ sinf cosdsingsin§ + cosBsin?6 + sinBsindcosdcosp. (10) 


In the same way we can find derivatives with respect to s in other directions, which 
are of interest in the deflection problem, 
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0 
a2 03 0.4 0.5 0.6 0.7 gauss 


Fig. 2. Deflection of w-mesons arriving at sea level from the zenith direction. 


dx/ds = —sinBsing cos6 + (cosfcosé — sinfsind cos) sin# (11) 
Knowing dx/ds and dh/ds, dy/ds is given by the relation 
(da/ds)? + (dy/ds)? + (dh/ds)? =1. 
Explicitly 
dy/ds = (sinBsin?6 cosy — cosf sind cosd) cos 6+ 
+ sinfsindsingsin§ + sinfcos?dcosp + cosfsindcosd. (12) 
We are also interested in the relative difference between the path length in the 


magnetic field and the path length without magnetic field. 
This is: 


eR BI _ int —tep sind 0086 os p- 


zt (cos 6 — tg sin 6 cos 6 cos n) | e0s 6ds+ tgP cosd sin p | sin 6 as| - (18) 
8 
4 6 


We may be interested in performing the integration: 


22 > 
. +2 
~ [e088 dasi— mss oats ° {cos Ods (12) 
2a) scosfB $ 0 
0 


where the dependence of f is included in cos. 
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s:cosB-h | s=30000m 
s-cos B W,= 0.1 GeV 
p=0 
9% 
5.05 


s=20000 m 
W,=0.1 GeV 


B=0 


s=20000m 
_-W1.0 GeV 

——"~  p=30° 

eae p= 270° 


s=20000m 

=1.0 GeV 
ae 
p=0",180 
s=20000m 
W,=1.0 GeV 
B=0° 


me $=20000m 
~~ W,=1.0 GeV 

B=30° 

=90° 


Fig. 3. Relative change of path length for a u-meson caused by deflection in a magnetic field. 


In fig. 2 the deflection of u~-mesons which arrive at sea level from the zenith 
direction and have started from a height of 20 km have been plotted as a function 
of 6 and B. The particle energy W, is 1 GeV at sea level. The values have been 
obtained by taking Af’ = arccos(dh/ds) from eq. (10), with 

f=), 

1 =7970m, 

0) = 0.001293 g/cm, 

q = 2(MeV cm?/g), 

W,= 90.107 GeV in formula (9) for 0. 


I 


The curves approximate very closely to straight lines. 

(scosB —h)/(scosB) have been plotted from eq. (13) as a function of B in fig. 3. 
The parameters are s, W,, 6 and ~. The dip angle 6 = 0. The values given in fig. 3 
justify the approximations made earlier by putting s = hcosf. 
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